AD  747772 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded 
by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person 
or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 


; :  _ _ _ 

!|.  13  tiVij  sr.:;:n  QV 

rlc  £.»  □ 


BV  . . 

ClSTBl-UT10H/AVAIt*SlLtTY  C03EJ 

! — Diii: — 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


\ 


* 


UNCLASSIFIED 


Sfomly  CUtsifieation 


DOCUMENT  CONTROL  DATA  •  R  &  D 

{Security  cjat  !Jic*tior^ot^lUOjJ*Hi/^>f^*b$jrMct^*n<Mnd*Kingjinnotrt^^ 


|f origin  a  ting  actiwitv  (Cotyormta  AufhorJ 

Naval  Postgraduate  School 
Monterey,  CA  93940 


2*.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


26.  GROUP 


N/A 


1 3  REPORT  TITLE 


Quasi  Area  Rule  for  Heat  Addition  in  Transonic  and  Supersonic  Flight  Regimes 


I  A  descriptive  NOTES  (Type  of  report  andjnclutiv  dataa) 

Final.  October  1970  to  August  19^1 


|  9  authoR(S)  (Flr*t  nmme,  middle  IniiiiT,  laat  name) 

Dr.  Allen  E.  Fuhs 


|0  REPORT  CATE 

August  19 72 


fli.  CONTRACT  OR  GRANT  NO. 
6.  PROJECT  NO.  ^066 

'•Task  306603 

d. 


KIFR  AP0-71-0ff> 


?*.  TOTAL  MO.  Or  PACE* 


129 


lb.  no.  or  per* 


24 


M.  O  PIC  IN  A  TOP'S  PEPOPT  NUMCCPIt) 

None 


•6.  OTHER  REPORT  NO(3)  (Any  othor  numkota  that  may  ba  aaolgnad 
thla  report) 

AFAPL-TR-72-10 


1*.  SUPPLEMENTARY  NOTES 


II.  A6STRACT 


1 0  OISTRIOU^ION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


112.  SPONSORING  MILITARY  ACTIVITY 


Air  Force  Aero  Propulsion  Laboratory 
Wright-Patterson  AFB,  Ohio  45433 


Body  shapes,  including  axisymraetric  and  three  dimensional,  have  been  developed  to 
minimize  wave  drag.  The  von  Karman  ogive  and  the  area  rule  are  examples. 

Similar  work  has  not  been  accomplished  for  optimum  shapes  with  propulsion. 

Propulsion  can  be  divided  into  two  categories — those  devices  with  internal  heat 
addition  and  those  with  external  burning.  For  internal  heat  addition  an  analytical 
model  is  formulated  which  introduces  the  propulsive  disc.  Attention  is  shifted  to 
external  burning,  which  is  eramined  for  one  dimensional  and  two  dimensional 
linearized  flow.  Heat  fronts  and  combustion  fans  are  discussed  as  examples.  Forces 
on  a  heat  source  in  a  uniform  stream  and  adjacent  to  bodies  are  derived.  Several 
possible  applications  are  examined  including  base  pressure  augmentation  by  external 
burning,  spin  recovery  using  external  burning,  and  transonic  boattail  drag 
alleviation.  Previous  work  on  base  pressure  augmentation  has  used  a  two 
dimensional  planar  model.  A  two  dimensional  axisymreetric  model  is  examined. 


J. 


DD .1473  IMGE  " 

S/N  0101-807-681* 


UNCLASSIFIED 

Security  CI***l(lCA«on 


»-SW0* 


. TrfliilT  • ' 


< 

AFAPLTR-72-10 


QUASI  AREA  RULE 
FOR  HEAT  ADDITION  IN  TRANSONIC 
AND  SUPERSONIC  FLIGHT  REGIMES 


Dr.  Allen  E.  Fuhs 
Professor  of  Aeronautics 
Naval  Postgraduate  School 


Approved  for  public  release;  distribution  unlimited. 


FOREWORD 


This  report  describes  work  accomplished  in  the  program  "Quasi 
Area  Rule  for  Meat  Addition  in  Transonic  and  Supersonic  Flight  Regimes" 
conducted  under  USAF  MIPR  APO-71-007.  The  work  was  accomplished 
during  the  period  1  Oct  1970  through  31  Aug  1971  at  the  Naval  Post¬ 
graduate  School,  Monterey,  California,  under  the  direction  of  Dr.  Allen 
E.  Fuhs.  The  report  was  submitted  on  1  Nov  1971. 

The  program  was  sponsored  by  the  Air  Force  Aero  Propulsion  Labora¬ 
tory,  Wright-Patterson  Air  Force  Base,  Ohio  under  Project  3066,  Tur¬ 
bine  Engine  Propulsion,  Task  306603,  Advanced  Component  Research. 

Dr.  Kervyn  D.  Mach,  AFAPL/TBC,  Turbine  Engine  Components  Branch, 
was  the  project  engineer. 

Publication  of  this  report  does  not  constitute  Air  Force  ap¬ 
proval  of  the  report’s  findings  or  conclusions.  It  is  published 
only  for  the  exchange  and  stimulation  of  ideas. 


ERNEST  C.  SIMPSON ^ 

Director,  Turbine  Engine  Division 


Body  shapes,  including  axisymmetric  and  three  dimensional,  have  been 
developed  to  minimize  wave  drag.  The  von  Karroan  ogive  and  the  area  rule  are 
examples.  Similar  work  has  not  been  accomplished  for  optimum  shapes  with 
propulsion.  Propulsion  can  be  divided  into  two  categories-- those  devices  with 
internal  heat  addition  and  those  with  external  burning.  For  internal  heat 
addition  an  analytical  model  is  formulated  which  introduces  the  propulsive  disc. 
Attention  is  shifted  to  external  burning,  which  is  examined  for  one  dimensional 
and  two  dimensional  linearized  flow.  Heat  fronts  and  combustion  fans  are 
discussed  as  examples.  Forces  on  a  heat  source  in  a  uniform  stream  and 
adjacent  to  bodies  are  derived.  Several  possible  applications  are  examined 
including  base  pressure  augmentation  by  external  burning,  spin  recovery  using 
external  burning,  and  transonic  boattail  drag  alleviation.  Previous  work  on 
base  pressure  augmentation  has  used  a  two  dimensional  planar  model.  A  two 
dimensional  axisymmetric  model  is  examined. 
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QUASI  AREA  RULE 

FOR  HEAT  ADDITION  IN  TRANSONIC 
AND  SUPERSONIC  FLIGHT  REGIMES 

Dr.  Allen  E.  Fuhs 
Professor  of  Aeronautics 
Naval  Postgraduate  School 

INTRODUCTION 

In  subsonic  flow,  drag  can  be  attributed  to  viscous  forces  on  the  body 
surface  and  to  the  release  of  vorticity  in  the  wake.  The  vorticity  is  associated 
with  the  generation  of  lift.  In  supersonic  flow,  there  is  a  new  phenomenon  lead¬ 
ing  to  drag;  it  is  the  radiation  of  waves  from  the  body.  These  waves  carry 
energy  whose  origin  is  the  work  done  by  drag  and  vehicle  motion. 

For  almost  every  case  of  a  body  flying  through  a  medium,  drag  is  considered 
to  be  a  degradation  of  performance.  Consequently,  considerable  effort  has  been 
devoted  to  optimum  shapes  giving  minimum  drag. 

One  approach  for  optimizing  a  body  shape  i3  to  represent  the  body  by  sources 
and  sinks. ^  For  a  closed  axisymmetric  body,  all  sources  and  sinks  are  on  the 
axis.  It  is  found  that  for  a  nonlifting  body,  the  slope  of  the  cross  sectional 
area  curve  equals  the  distribution  of  sources  and  sinks,  f(x).  Distance  along 
the  axis  is  x.  By  means  of  the  calculus  of  variations,  one  can  arrive  at  shapes 
giving  minimum  wave  drag.  There  are  constraints  on  the  body  geometry,  and  the 
optimum  shape  depends  on  these  constraints.  Two  examples  are  given  in  Table  I. 
The  shapes  given  in  Table  I  are  for  minimum  wave  drag.  Viscous  effects  are  not 
included,  and  the  bodies  are  nonlifting. 


Table  I.  Optimum  Shapes 


Name  of  Body 

Constraints 

Wave  Drag 

( 2 ) 

Von  Karman*  Ogive'  ' 

Given  length,  l ,  and  base  area,  S 

4  /  sf\ 

M  l2) 

Sears**Ogive^ 

Given  length,  l ,  and  volume,  V 

128  / 

; q|  t*j 

•Used  on  Low  Altitude  Short  Range  Missile. 
••Discovered  by  W.  R.  Sears  of  Cornell  University. 


It  was  discovered  by  Lomax  and  Heaslet^  and  by  Whitcomb^  that  the 
optimum  3hape  for  an  axisymmetric  body  applies  to  slender  bodies  with  wings  and 
rudders  if  one  interprets  the  area  distribution  properly.  It  is  the  area  cut  by 
the  Mach  cone.  If  one  uses  this  area  distribution  from  a  three-dimensional  body 
to  match  the  area  distribution  of  an  optimum  axisymmetric  body,  then  minimum  wave 
drag  occurs.  Thi3  procedure  is  known  as  the  area  rule. 

For  a  body  with  thrust,  there  are  several  nev  variables  introduced.  Thrust 
may  be  developed  by  adding  heat  to  a  flow  internal  to  the  body.  Thrust  also  may 
be  obtained  by  adding  heat  external  to  the  body. 

Consider  external  burning.  Is  there  an  area  rule  for  a  body  with  external 
heat  addition  analogous  to  the  area  rule  involving  body  shape  onlyT  For  this  case, 
there  are  two  functions  to  optimize.  One  is  body  shape,  and  the  other  is  the 
distribution  of  heat  addition.  Table  I  gives  some  typical  constraints;  an 
additional  constraint  for  the  optimization  could  be  specific  fuel  consumption. 
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Instead  of  finding  minimum  drag,  the  procedure  should  give  the  best  thrust- 
minus-drag  or  zero  thrust-minus -drag  depending  on  the  constraints.  Since  both  the 
distribution  of  heat  and  cross  sectional  area  are  involved,  an  appropriate  name 
is  "quasi  area  rule." 

When  thrust  is  developed  by  adding  heat  to  a  flow  internal  to  the  body,  the 
quasi  area  rule  suggested  in  the  previous  paragraphs  does  not  apply.  For  in¬ 
ternal  burning,  the  waves  due  to  heat  addition  obviously  do  not  occur  as  in  the 
case  of  external  burning.  There  is,  however,  a  change  in  stagnation  pressure 
and  temperature  between  the  inlet  and  nozzle,  Also  the  velocity  vector  changes 
from  inlet  to  exhaust.  An  energy  disc  can  represent  those  changes.  This 
report  discusses  aspects  of  the  thrust-minus -drag  problem  for  an  airbreathing 
engine  with  confined  heat  addition.  Figure  1  illustrates  the  various  problems 
that  car.  be  studied  with  internal  burning  when  the  emphasis  is  on  thrust-minus- 
drag.  The  analysis  is  linearized.  Two-dimensional  (2D)  planar  is  the  geometry 
of  a  propulsive  wing.  Axisymmetric  2D  geometry  is  an  isolated  nacelle  without 
fins  or  wings.  Addition  of  fins  or  wings  yields  the  2D  propulsion  with  3D  body. 

The  connection  between  these  two  geometries  is  the  area  rule. 

Figure  2  outlines  the  problems  of  interest  for  external  burning.  The 
analysis  is  linearized.  For  the  2D  planar  case,  the  heat  and  mass  sources  can 
be  distributed  over  a  surface  or  throughout  a  volume.  Volume  distributions  must 
meet  requirements  for  slender  body  theory.  Axisymmetric  2D  bodies  are  related 
to  3D  bodies  by  the  quasi  area  rule. 

The  linearized  2D  planar  flow  with  heat  addition  has  been  studied  extensively. 
A  solution  to  the  2D  axisymmetric  case  does  not  appear  in  the  scientific 
literature.  Solutions  for  3D  bodies  with  heat  addition  have  not,  as  yet,  been 
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Figure  2.  Classification  of  Flows  with  External  Burning 
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obtained.  In  order  to  develop  the  quasi  area  rule,  these  solutions  are 
essential  building  blocks. 


PROPULSION  WITH  INTERNAL  HEAT  ADDITION 


Introduction 


For  a  supersonic  aircraft,  the  variation  in  nozzle  pressure  ratio  (NFR) 
is  very  large.  This  requires  a  variable  geometry  nozzle.  Under  certain  flight 
conditions,  the  nozzle  setting  and  nozzle  flow  result  in  large  base  or  boattail 


drag. 

One  motivation  for  treating  the  propulsion-airframe  problem  as  a  unit  from 
the  area  rule  point  of  view  is  to  obtain  new  design  approaches.  First,  the 
typical  variation  in  NPR  and  nozzle  geometry  will  be  calculated  for  a 
Mach  2.5  aircraft.  These  results  will  be  discussed  to  illustrate  the  problem. 
Second,  the  model  for  representation  of  a  propulsion  system  with  internal  heat 
addition  will  be  developed.  Third,  approaches  to  solution  of  the  model  will  be 
outlined. 

Exhaust  Variation  with  Mach  Number 

The  aim  of  this  section  is  to  calculate  nozzle  pressure  ratio,  nozzle  area 
ratio,  mass  flow  ratio,  and  thrust  of  an  ideal  turbojet  as  a  function  of  Mach 
number.  Knowledge  of  the  trends  is  essential  in  understanding  aircraft  exhaust 
problems. 

Consider  an  ideal  turbojet  with  an  ideal  inlet.  There  is,  however,  a 
restriction  on  the  Mach  number,  Mg,  into  the  compressor.  The  thrust  is  given  by 

I  F  \/&og\ 


F  \V/wdg  (1) 

where  ma  is  a  nu  flow  of  air,  and  fhj  is  mass  flow  of  inlet  at  design  conditions. 
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For  the  inlet,  the  various  areas  are  illustrated  in  Figure  3(a).  Station  0 
is  at  a  point  in  the  freestreara  where  i3  ambient  pressure.  Stations  1  and  .? 
are  inlet  throat  and  compressor  face  respectively.  The  mass  flow  ratio  is 


m 


M0  A0  At 


M, 


Od  *t 


At  Al 


(2) 


where  Mq  is  flight  Mach  number,  and  M^  is  design  flight  Mach  number.  Areas 


Od 


Aq<  A^»  and  Aj  are  identified  in  Figure  3(a).  Both  a  fixed  area  and  variable 
geometry  inlet  are  considered.  The  ratio  Aq/A^  is  obtained  from  isentropic 
flow  tables  for  given  Mq  for  the  variable  geometry  inlet.  For  a  fixed  geometry 
inlet  Aq/A^  comes  from  flow  tables,  entering  the  flow  tables  with  the  Mach 
number  downstream  of  a  normal  shock  at  Mq.  The  mass  flow  ratios  are  shown  in 
F igure  4  for  both  cases.  A  variable  geometry  operates  supercritically  and 
swallows  more  air.  A  variable  geometry  ramp  or  spike  inlet  may  have  superior 
mass  handling  capability  at  subsonic  speeds. 

The  thrust  per  pound  of  air  for  an  ideal  turbojet  i3  given  by 


f  !o”o  n~ 

*  V7^"1 


- 1 


(3) 


where  the  symbols  not  previously  defined  have  the  following  meaning: 

Sq  *  freestream  speed  of  sound 

Tr  =  T«pq/Tq  *  ratio  of  stagnation  to  static  freestream  temperature 
Tc  -  =  3^a6na^'i°n  temperature  ratio  across  the  compressor 

The  symbol  T*  deserves  more  discussion.  It  is 


T*  =  Wb 


TT4^T0 


(4) 


where  T ^  =  ^T4^T3  =  stagnation  temperature  ratio  across  the  burner  or  combustor. 

The  value  of  T^  is  limited  and  has  a  maximum  value.  Hence  for  a  given  engine 

and  altitude,  T*  is  fixed  when  the  engine  is  at  maximum  power  or  maximum  TIT. 

7 


®  ©  ® 


A*  LIP  AREA 

At  THROAT  AREA 

Ac  COMPRESSOR  FACE  INLET 


(a)  AREAS  FOR  INLET 


I  till  III 
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(b)  STATION  NUMBERING 


FIGURE  3.  NOTATION  FOR  TURBOJET  CALCULATIONS. 
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The  station  numbering  is  shown  in  Figure  ^(b).  The  throat  area  for  the 


nozzle  is 


gr sA 

r  ptr 


neglecting  fuel  addition.  R  is  the  gas  constant.  P  is  a  function  of  the  ratio  of 
heat  capacities,  y,  and  has  a  value  of  .683  for  y  =  1.4.  The  nozzle  pressure  ratio, 
NPR,  is  p^,g/p^.  A  ratio  of  ky  to  the  design  value  is 

A~  ra  NFR.  /t_- 

—  ■  j  — S-.pL  (6) 

A7d  "d  NPR  V  TT8d 

The  area  Ag  is  obtained  from  NPR,  isentropic  flow  tables,  and  the  assumption 

pg  =  Pq.  Two  useful  ratios  are  Ag/Ac  and  A^/Ac,  i.e.,  comparing  nozzle  areas  with 

compressor  face  area. 

Calculations  were  made  for  the  values  shown  in  Table  II. 

One  of  the  results  of  the  calculation  is  NFR  which  is  shown  in  Figure  5.  Some 
of  the  other  results  are  summarized  in  Table  III  The  symbols  for  the  turbine  ratios 
are  \  =  T^/T^  and  *t  =  PT5/pT4. 

Figure  6  illustrates  the  area  variations  required  for  a  supersonic  turbojet 
operating  at  maximum  TIT  for  different  Kq.  Examination  of  Figure  6  and  the  data  in 
Table  III  shows  why  turbojet  nozzles  usually  operate  underexpanded  at  design  Mach 
number.  The  percentage  variation  in  A7  is  38  per  cent  based  on  design  A^.  The 
percentage  variation  in  Ag  is  77  per  cent  based  on  design  Ag.  NPR  varies  from  4  at 
low  subsonic  flight  to  more  than  25  at  design  0f  2.5. 
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Table  II.  Values  for  Calculation  of  Turbojet  Quantities 

Maximum  turbine  inlet  temperature 

TT4 

2000°F 

Altitude 

Sea  Level 

Freestream  density 

Po 

.002378  slugs/ft3 

Freestreara  speed  of  sound 

a0 

1117  ft/sec 

Compressor  area 

Ac 

8.05  ft2 

Tip  radius  at  face 

r* 

O 

2  ft 

Hub  radius  at  face 

rh 

1.2  ft 

Maximum  Mach  number  at  compressor 

M 

n2roax 

0.65 

Inlet  throat  area  at,  design 

Atd 

7.1  ft2 

Design  flight  Mach  number 

M0d 

2.5 

Compressor  pressure  ratio  at 

100  per  cent  corrected  speed 

n 

c 

8 

tfSess 


EXIT  OF  NOZZLE  FOR 


As  shown  in  Figure  6,  to  obtain  ideal  expansion  the  nozzle  is  converging- 
diverging  for  all  Mq  illustrated.  Usually  the  nozzle  is  a  simple  converging 
nozzle  from  subsonic  through  transonic  range  of  Mq.  This  results  in  a  large 
boattail  area.  At  the  design  MQ  the  exit  area  for  ideal  expansion  significantly 
exceeds  the  maximum  area  of  the  propulsion  system.  The  nozzle  is  cut  off,  relative 
to  ideal,  at  some  area  appropriate  to  maximum  engine  area.  The  result  is  under¬ 
expansion  and  a  ballooning  plume.  A  shock  wave  may  originate  at  the  diverging 
streamlines  of  the  plume,  giving  a  drag  increment. 

Model  for  Propulsion  with  Internal  Heat  Addition 

Consider  a  propulsion  device  as  sketched  in  the  top  of  Figure  7.  To 
represent  the  wave  drag  of  the  inlet,  it  is  necessary  to  duplicate  streamline 
A-A.  It  does  not  matter  what  the  details  are  of  the  inlet  geometry  internal 
to  the  streamtube  bounded  by  A-A. 

The  inner  body  which  generates  a  flow  having  streamline  A-A  is  shown.  Note 
that  it  is  not  the  same  shape  as  the  inlet  spike.  A  distribution  of  sources  on 
the  axis  could  generate  the  inner  body  as  well  as  streamline  A-A.  Note  that  the 
spike  shock  wave  does  not  have  the  same  angle  as  the  cowl  shock  wave. 

If  one  applies  the  momentum  relations  to  the  control  volume  of  Figure  7, 
the  thrust  minus  drag  is  given  by 


15 


FIGURE  7.  MODEL  OF  AN  INLET  USING  SOURCES. 


where  u^»  and  u^.  are  the  dimensionless  perturbation  velocities  arising 
from  the  nacelle.  It  is  apparent  that  for  the  case  where  the  exhaust  exit 
defines  a  plane,  as  illustrated  in  Figure  7,  the  thrust  and  drag  can  be 
clearly  separated.  Waves  generated  by  the  exhaust  plume  do  not  influence  T  -  D. 

Following  the  analysis  suggested  by  the  model  of  Figure  7,  one  can  obtain 
a  nacelle  of  known  wave  drag.  Figure  8  illustrates  the  concept.  Start  with  an 
ogive  of  known  shape.  Plot  the  streamlines  in  the  flow  adjacent  to  the  ogive. 
Pick  the  radii  r^  and  Tg  suitable  for  the  nacelle.  Insert  the  nacelle  within 
the  streamtube.  The  wave  drag  of  the  ogive  and  nacelle  will  be  the  same.  The 
drag  can  be  calculated  using  the  potential  function  for  the  ogive.  In  the 
lower  part  of  the  figure,  a  nacelle  with  a  spike  inlet  is  shown. 

Consider  some  thrust  device  inside  a  momentum  control  volume  shown  in 
Figure  9.  How  does  one  modify  the  variables  so  as  to  change  conditions  on  the 
surface  of  control  volume  and  improve  T  -  D?  What  doe3  one  have  available  to 
change  the  wave  drag?  Available  variables  are:  (a)  number  and  location  of 
nosales,  (b)  velocity  profile  at  nozzle  exit,  and  (c)  body  geometry. 

To  obtain  an  optimum  conf iguration,  one  must  have  a  method  to  represent 
the  body  and  exhaust  analytically.  The  analytical  representation  can  be 
subjected  to  the  calculus  of  variations. 

The  body  can  be  represented  by  sources  and  sinks.  Changes  in  stagnation 
pressure,  stagnation  temperature,  velocity  vector,  and  other  flow  quantities 
make  representation  of  exhaust  more  difficult. 

Carrying  the  idea  of  Figure  7  further,  let  us  look  at  an  ideal  ramjet. 

For  an  ideal  ramjet,  the  following  equalities  hold  true: 


) 

*  FIGURE  9.  MOMENTUM  CONTROL  VOLUME  FOR  A  THRUSTING 

BODY. 


19 


"“Jmi  ItatsJi  A  ’  Z  i  li  XKiA&s ,1 M ! 


*1  =  A6 


PT6  =  PT1 


v6  _  TT6 
vo  TT1 


p6  _  TT1 
P1  TT6 


pl  =  p6 


For  an  ideal  ramjet  (or  for  any  other  propulsion  device  satisfying  the  above  set 
of  equalities),  one  can  represent  the  flow  by  sources,  sinks,  and  an  "energy 
disc."  This  is  illustrated  in  Figure  10;  stations  1  and  6  are  identified  in 
Figure  10.  The  inner  body  duplicates  streamline  ABC.  The  inner  body  is  formed 
from  the  source-sink  distribution  in  the  right-hand  column  of  Figure  10.  The 
change  in  density,  velocity,  and  stagnation  temperature  is  achieved  by  an 
"energy  disc."  In  fact,  that  defines  the  energy  disc.  It  changes  the  flow  so  that 
the  equalities  list  above  is  correct.  Tne  sources,  sinks,  and  energy  disc  shown 
in  the  bottom  of  the  right  column  would  have  the  same  thrust  and  drag  as  the 
ideal  ramjet  sketched  in  the  left-hand  column. 

This  model  looks  a3  if  it  is  satisfactory  for  an  ideal  ramjet,  but  it  has 
some  drawbacks; 

1.  The  energy  disc  is  not  a  heat  front  or  other  simple  fluid  flow. 

2.  When  /  A^,  it  may  be  difficult  to  apply  the  energy  disc. 

3.  For  the  case  p^  4  P^»  some  means  of  changing  p^,  must  be  added. 

The  energy  disc  is  bounded  by  a  streamtube.  Changing  p,p  within  the  disc  gives  a 
mismatch  across  the  disc  in  regard  to  streamline  shape. 

A  satisfactory  way  of  specifying  the  conditions  at  station  6  is  needed,  and 
the  energy  disc  has  been  suggested.  One  wants  to  retain  sources  and  sinks  for 
two  reasons.  Wave  drag  is  readily  calculated  if  you  know  the  distribution  function. 
The  distribution  function  can  be  optimized  to  yield  minimum  wave  drag. 
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FIGURE  K>.  IDEAL  RAMJET  REPRESENTED  BY  SOURCES,  SINKS,  AND  AN  ENERGY  DISC. 


The  energy  disc  represents  the  propulsion  device.  A  particular  engine 
generates  changes  in  p,p,  T^,,  and  velocity  between  the  inlet  ana  exhaust.  For 
the  case  of  a  plane  for  an  exit  as  shown  at  station  6  in  Figure  10,  one  does  not 
need  the  complexity  of  an  energy  disc.  Equation  (7)  indicates  that  drag  and 
thrust  are  separable  for  this  case.  When  one  varies  quantities  (a),  (b),  and 
(c)  discussed  previously  as  available,  independent  variables,  then  the  energy 
disc  may  be  an  attractive  model. 

This  section  is  closed  by  noting  that  a  model  has  been  proposed  but 
has  not  been  developed, 

EXTERNAL  BURNING;  BASIC  EQUATIONS  AND  SOLUTIONS 
One-Dimensional  Heat  Addition 

To  provide  a  basis  for  comparison  with  subsequent  sections,  one-dimensional 
heat  addition  will  be  briefly  reviewed. 

There  are  several  approaches  to  one-dimensional  heat  addition.  For  unsteady 

flow  there  is  the  3elf  similar  technique  developed  by  Sedov. ^  ye  consider 

steady  flow  only  here.  There  is  the  famous  solution  for  constant  area  heat 

addition,  i. e. ,  the  solution  by  Rayleigh.  Rayleigh  flow  is  discussed  in  most 

(7) 

gas  dynamics  books.  When  there  is  simultaneous  variation  of  area  and  heat 
addition,  the  techniques  outlined  in  the  Princeton  Series,  Volume  III,  are 
applicable. ^ 

(9) 

The  influence  coefficients  developed  by  Shapiro  and  Hawthorne  will  be 
discussed  here.  Table  IV  summarizes  the  influence  coefficients  for  heat 
addition  and  isentropic  flow  with  area  change.  The  influence  functions  relate  a 
fractional  change  of  a  dependent  variable,  say  M2,  to  a  fractional  change  of  the 
independent  variable,  A  or  h^,.  So  far  as  heat  addition  is  concerned,  the 


22 


^’^*^®SSW5SB@»' *tS\.rHt.~n***4s£*  *8&&$&&Ei 

l^^Ws«*r¥maw*^  **~,*^.t  _ _ .... 

31 


f?l.^1,<'--5ft-4^’*W-*'*d<  ■*'*- 


Table  IV. 

Influence  Coefficients  for  Area  Change  and  Heat  Addition 

dA 

dh 

T 

Tf 

P 

dM2 

-2(1  ♦  ^-—M2) 

1  +  yM2 

M2 

l-K2 

1  -  M2 

& 

1 

1 

u 

1  -  M2 

1-M2 

dT 

(v  -  Dm2 

1  -  YM2 

T 

1  -M2 

1-M2 

d£ 

M2 

1 

P 

1  -  M2 

i 

(— i 

yM2 

P 

1  -  M2 

1-M2 

dF 

1 

n 

F 

i  +ym2 

ds 

c 

0 

i 

P 

velocity,  u,  and  static  pressure,  p,  change  as  they  would  flow  into  decreas¬ 
ing  area.  For  u  and  p,  heat  addition  and  area  variation  produce  the  same 
result. 

For  unconfined  flow  with  heat  addition,  there  are  some  analogies  between 
the  flow  due  to  a  particular  body  shape  and  the  flow  due  to  a  given  distribution 
of  heat  addition.  These  analogies,  which  are  limited  in  scope,  will  be  pointed 
out  as  they  arise  in  subsequent  discussion.  The  fact  that  heat  addition  and  a 
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body  may  influence  the  flow  the  same  way  may  be  understood  using  the  du/u  and  dp/p 
influence  coefficients. 

The  influence  coefficients  are  the  first  terms  in  a  Taylor  series  expansion 
of  the  dependent  variables  in  terms  of  the  independent  variables.  Many  of  the 
results  of  linearized  heat  addition  in  an  unconfined  flow  can  be  related  to 
influence  coefficients.  This  statement  is  applicable,  for  examble,  to  heat 
fronts  and  combustion  fans. 

Linearized  Equations  for  Combined  Mass  and  Heat  Addition 

The  equations  of  motion  for  steady  flow  with  heat  and  mass  addition  are  as 
follows  J 


Continuity 

Momentum 


div(pu )  =  m 


Du  -*■  . 

pDt  +  VP  +  =  0 


(8) 

(9) 


Energy 


div  piT(h  +  4”)  =  Q 


(10) 


In  the  preceding  equations,  m  is  the  rate  at  which  mass  is  added  per  unit 
volume.  Likewise  Q  is  the  rate  at  which  heat  is  added  per  unit  volume. 

To  linearize  these  equations,  one  introduces  perturbation  quantities  vhleh 
are  defined  as  follows: 


P  =  +  P'  ;  P  =  P«,  +  P‘ 

u  =  e  (U  +u')+^V  ♦  e  w ’ 
x  y  z 


+  T»* 


(11) 


When  the  perturbation  quantities  are  introduced,  the  resulting  equations  are 


Continuity  p  u’  ♦  p  +  Up^.  =  m 

##  X  ir  v 


Momentum 


r<*  y  'x 


P„u«ux  +  Px  *  "  mU® 

P*V»  *  Py  =  0 


(12) 

(11) 


24 


Energy 


P«P«  a/p' 

r  -  1  8x\P#0 


TP' 

71 


(—  -  "->  -  Q 


(14) 


We  have  formulated  a  2D  planar  flow  so  that  w'  =  0.  To  avoid  writing 
subscript oo,  ve  now  drop  it  in  subsequent  equations.  The  notation  with 
subscript  x  means  partial  differentiation  with  respect  to  x. 

Starting  with  Equations  (8)  to  (10),  one  could  obtain  the  linearized 
equations  appropriate  for  heat  and  mass  addition  in  axisymmetric  flow.  The 
equations  have  been  solved  for  mass  addition  (sources)  but  jt  for  heat 
addition  in- axisymmetric  flow. 

The  solution  for  a  source  in  two-dimensional  planar  flow  will  now  be 
obtained.  Consider  a  row  of  point  sources  evenly  distributed  along  the 
z-axis.  Axis  orientation  relative  to  the  flow  is  shown  in  Figure  11.  For 
a  point  F(x  -  Xp  y  -  y^  ),  only  the  sources  within  the  Mach  forecone  influence 
po  Lnt  P  as  shown  in  Figure  11.  To  find  the  potential  in  two-dimensional  flow, 
add  up  the  sources  along  the  z-axis  from  k  to  d.  Points  k  and  d  are 
illustrated  in  Figure  11.  From  the  geometry  of  Figure  11 

zd =  ~  xr  tan  m  “  (y  “ 

Integration  of  a  row  of  point  sources  yields  the  potential 

z. 


(15) 


0(x,y,O)  =  - 


«  «  • 

-  z . 


dz. 


2"  -y /(x  -x1)2  -  ?z\  -  e2(y  -  yjY 


(16) 


This  integrates  to 
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FIGURE  II.  POTENTIAL  FOR  A  LINE  SOURCE  IN  SUPERSONIC  FLOW. 
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0(x, z)  =  -  —  sin" 
fn 

Inserting  the  limits  gives 


'(x  -  y^)2  -  f2(y  -  yi)2 


0(x»y)  * 


The  potential  is  a  constant.  Think  of  a  series  cf  Mach  cones  with  apex 
at  the  z-axis.  Upstream  of  these  cones,  the  potential  is  zero.  Downstream 
of  the  envelope  of  Mach  cones,  it  is  constant.  One  can  represent  this  by 


0(x,y)  =  "  2$  I  ['*  *  *1^  -  ”  yl^] 

where  I  is  the  unit  function.  The  velocity  components  are,  for  a  unit 
source  strength, 

u  =  =  "  2gS[(x  “  xl)  *  p(y  “  yl|) 

and 

w  =  =  +  l6[(x  "  O  -  P(y  -  yx)] 

The  (♦)  sign  in  front  of  p  gives  right-running  waves,  whereas  the  (-)  sign 
gives  left-running  waves. 

The  solutions  for  a  heat  source  in  supersonic  flow  are^^ 
u’  *  -  6(x  ■ p3r) 

v  ■  -kjj^SOc  -  (7) 

?•  -  -  (r) 

p.  -  &  %A)Ha-  fiU  -  sjt)  -  S(y)l(x) 

2a  p  a  U 


One  can  verify  that  there  are  the  appropriate  solutions  when  m  *  0  and  Q  /  0 
by  substitution  into  Equations  (12)  to  (14).  The  second  term  in  Equation  (25) 
is  a  wake.  Pressure  in  the  wake  is  the  ambient  value,  and  the  velocity 
returns  to  freestream  value.  The  solutions  in  Equations  (22)  to  (25)  are 
valid  for  all  x  and  y  including  the  source  point  x  =  y  a  0.  The  symbol  q, 
which  has  units  heat/(length)(time),  will  be  discussed  shortly.  The  line  heat 
source  in  a  planar  flow  is  illustrated  in  Figure  12, 

Similar  solution?  can  be  written  for  the  case  of  mass  addition;  these  are: 

u’  =  “  |p  6(x  -  for)  (26) 

V=^6(x-py)  (2.7) 

p'  =  fjj3  6(x  -  0y )  (28) 

P’  B  ffS“$(x  -  for)  (29) 

Substitution  of  Equations  (26)  to  (29)  into  Equations  (12)  to  (14) 
indicate  these  are  solutions  except  at  the  mass  source  point  x  =  y  =  0. 

With  Equations  (22)  through  (29),  one  has  solution  for  volume  distribution 
of  heat  and  in  regions  where  ra  =  0.  This  is  adequate  to  describe  planar 
bodies  with  adjacent  heat  addition  zones. 

Volume  Heat  Addition 

Since  the  heat  may  be  added  throughout  a  volume  rather  than  a  line  source 
(the  z-axi3  in  Figure  11),  it  is  worthwhile  to  formulate  heat  a3  a  volume 
release.  Define,  as  in  Equations  (10)  and  (14) 


L span M area  Msec, 


Then  we  define 

heat 

^  ~  (span) (sec) 
q  =  Qdxdy 

Another  area  element  may  be  better  than  dxdy.  Some  element  with  its  sides 
parallel  to  a  Mach  line  is  preferred.  We  think  of  Q  =  Q(x,y).  Figure  13 
is  a  series  of  drawings  illustrating  the  variables. 

To  get  heat  per  unit  span 


H  = 


'  0  **  0 


Qdxdy  =  QU 

L 


(’0) 


if  Q  =  constant 


Heat  released  in  the  complete  volume  is  It  =  Hs. 

heat  „  heat 


U’ 


H  = 


sec  (span) (time) 

Now  to  tie  into  the  2D  flow  equations,  we  use  Equation  (22)  as  an  example, 

u*  *  -  6(x  "  fy) 

At  a  point  on  •  characteristic,  u'  results  from  a  source  on  that  same 
characteristic.  If  there  is  more  than  one  source  on  that  characteristic, 
one  adds  the  influence  of  each. 

(y  -  l)(q1  ♦  q2) 


(31) 


u 


I  c 


2ypP 


6(x  -  py) 


(32) 


For  many  heat  sources,  in  fact  n  heat  sources, 

i  *  n 


u,  .t  . 

2rp? 


5(x  -  py) 


(33) 
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FIGURE  13. 

ILLUSTRATION  OF  QUANTITIES  RELATED  TO  SPECIFICATION  OF  HE  AT  RELEASE. 


The  equation  for  u*  is  correct  if  one  thinks  of 


q*  3  Qdy  = 


l  span  nlength  H  sec , 


when  the  heat  source  is  distributed  in  the  volume.  When  heat  sources  are  not 
discrete  but  are  distributed,  then  one  wants  to  integrate  along  a  characteristic. 
Let  s  be  distance  along  the  characteristic.  We  express  in  terms  of  Q  and  ds 
»,  qi  =  Q  sin  ds  =  dq  (34) 

The  perturbation  velocity  is 


f  3  r  3 

u'  =  “  ^2^  I  Q  sin  ji  ds  =  -  Q(x,y)ds 

*  j  . 


Since  s,x,y  are  related  quite  simply 

x  -  s  cos  jx,  y  -  s  sin  jxt  x  =  py 

one  can  integrate  u'  several  ways: 

f3  f7 

sin  p  I  q|s  cos  jif  s  sinjujds  or  I  Q(fy,y)dy 


or  i  Q(x,x/p)dx/p 


Referring  to  the  region  bounded  by  0  <  x  <  L  and  0  <  y  <  l  in  Figure  13(b) 
and  for  case  Q  *  constant  inside  the  region 
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For  a  point  inside  the  heat  release  region 


X1  «  X, 


rWAKB'~I 


(xT,y)  =  - 


Q(x',y)dx' 


The  symbol  x^  is  defined  in  Figure  14. 

Heat  Addition  at  a  Plane:  A  Heat  Front 

Consider  a  distribution  of  line  sources  which  form  a  plane  at  angle  0 
to  the  flow,  as  illustrated  by  Figure  15.  Terminology  similar  to  that  as*-'.'’ 
in  supersonic  wing  theory,  i.e. ,  subsonic  and  supersonic  ''leading  edge,"  1* 
employed  to  identify  two  cases,  9  ^  ju  or  Mft  ^  1.  The  normal  Mach  number  Mq 
is  M  sin  0.  Since  sin  y  *  1/M,  the  normal  Mach  number  is 

V®?  <«> 

You  can  see  from  an  examination  of  Figure  15  that  the  heat  source  can  influence 

flow  upstream  if  M  <  1.  However,  if  M  >  1,  the  heat  source  does  NOT  influence 
n  n 


flow  upstream  contrary  to  statements  made  by  Tsien  and  Beilock. 


A3,  so  it 


is  apparent  that  point  P  is  influenced  only  by  right-running  characteristics  for 

M  <  1,  whereas  P  is  influenced  by  both  families  of  characteristics  for  M  >  1. 
n  n 

The  perturbation  velocities  for  a  heat  front  will  now  be  obtained. 

Using  Equation  (22)  and  integrating  along  the  heat  front,  the  formula  for  u*  is 

oo 

r 

u'(x,y)  -  -  -^—^3  g  [(x  -  s  cos  0)  -  $(y  -  s  sin  Q)|  ds  (46) 


where  s  is  distance  along  the  heat  front.  Consequently, 


~~ ■zzszmssi+z 


M 


p 


$ 


-  1 


u.  =  k-zJ-h  — 

2ypp  p  ain  0  -  cos  0 
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coa  0  -  p  sin  0 


;]da 


(47) 


- ,.....(t..- l-),q - - 

-VpPCP  sin  0  -  cos  0) 


Since  p  =  ^M2  -  1,  it  is  also  true  that  p  =  l/tan  ju.  Using  this  fact 

..»  „  (y  -  l)q  tan  n  sin  n 
2yp  sin(p  -  0) 


(48) 


For  left-running  waves,  p.  is  positive;  and  for  right-running  waves,  p.  is 
negative.  For  a  subsonic  "leading  edge,"  0  is  less  than  p.  Downstream  of  the 
heat  source  for  Mn  <  1,  Equation  (48)  indicates  u^  will  be  negative;  upstream 


of  the  source  u^  will  be  positive.  In  addition 

kl  >  kl 


(for  Mft<  1) 


where  subscripts  u  and  d  indicate  upstream  and  downstream  respectively. 

For  a  supersonic  "leading  edge,"  which  is  illustrated  on  the  right-hand 
side  of  Figure  15,  the  left  running  Mach  wave  gives  a  u£  which  is  negative. 

The  right-running  Mach  wave  gives  a  u^  which  i3  also  negative  but  which  is  less 
in  absolute  value  than  u£.  Subscripts  L  and  R  denote  left-  and  right-running 
waves.  Consequently, 

“d  ’  “i  *  “h  <0 

Now  let's  look  at  v';  from  Equation  (23)  one  obtains  for  this  case* 

^  co 

6  |(x  -  s  cos  0)  -  p(y  -  s  sin  9)J  ds 

'-CO 


v«(x  y)  =  ir. r  Ik 

V  KX,/J  2yp 


(49) 
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Integration  of  the  6 -function  yields 


vi  =  (Y-r-Uq _ 1 _ _  (L-r-Uq - 4n-E_ 

2yp  P  sin  9  -  cos  0  2yp  sin  (0  -  p) 

For  the  case  Mn  <  1,  is  negative,  and  v^  is  also  negative. 

Furthermore 


(50) 


tx  s 
u 


(for  Mn<  1) 


For  the  case  Mq  >  1,  v£  is  positive  while  v^  is  negative.  In  add:  v.ion 

vd  3  VL  ‘  ki  >0  (for  Mn>  1} 

which  means  the  flow  is  turned  toward  the  line  source  much  the  same  as  an 
oblique  shock  wave.  For  Mn  <  1,  the  flow  turns  away  from  the  heat  front. 

A  natural  question  i3,  ’’Is  there  any  special  direction  to  the 
perturbation  velocity  vector?"  Let  “v  be  the  resultant  perturbation  velocity 
vector.  It  is  given  by 


v 

r 


tan  u  sin  u 
sin(©  -  ji) 


+ 


tan  n  sin  u\  •+  r  sin  u  _  sin  u 
sin(0  ♦  pi)  y  \iin(9  -  p)  “  sin(9  ♦  p) 


(51) 


The  algebraic  signs  are  correct  for  the  case  Mfi>  1  which  means  9  >  p.  The 
symbol  C  equals  (y  -  l)q/2yp.  Let  a  be  the  angle  of  the  perturbation 
velocity  vector  relative  to  freestream  direction.  Then 


sin  u 


sin(0  -  n) 
tan  a  =  - 5 - 


sin  n 
sin (9  ♦  u) 


tan  u  sin  »  ,  tan  u  sin 
sin(©  -  n)  sin(©  '*  fi) 


Manipulation  of  Equation  (52)  leads  to 


tan  o 


_1 _  sin  u  cos  9  _ _ 1 


tan  fx  cos  p  sin  9  tan  9 


(52) 


(51) 


The  conclusion  from  Equation  (53)  is  that  the  velocity  vector  is  normal  to 
the  planar  heat  source.  Since  this  i3  so,  one  can  transform  the  supersonic 
case  of  heat  addition  to  an  equivalent  ID  problem.  To  do  thi3,  one  uses 


(54) 


du  (xa£  ud)2  ♦  (v^  +  vp2 

— -  JC  — . . . . — —  — »»■" 

u  u 

where  u  =  U  sin  8.  Proceeding  by  substitution  of  components  i^om 
Equation  (51)  into  Equation  (54)  gives 


du  1  hT 

U  sin  C  1  -  M2  c  T 
n  p 


(55) 


which  is  precisely  the  influence  coefficient  result  when  properly  interpreted. 
See  Table  IV.  For  a  subsonic  "leading  edge,"  the  velocity  increases.  This 
means  that  the  flow  is  turned  away  from  the  planar  heat  source  for  Mr  <  1 
(analogous  to  Prandtl  Meyer  expansion)  and  toward  the  planar  heat  source  for 
Mn>  1  (analogous  to  oblique  shock  wave). 

As  has  been  demonstrated,  one  of  the  ways  of  adding  heat  which  lends 
itself  to  analysis  is  by  means  of  a  heat  front.  The  example  chosen  by 
Tsein  and  Beilock^^  was  a  heat  front.  Oswatitsch^  ^  considers  other 
aspects  of  heat  fronts. 

Heat  addition  at  a  front  is  comparable  in  complexity  to  linearized  heat 

addition.  Combustion  is  a  volume  phenomenon.  However,  if  the  length  of  the 

heat  release  region  is  much  less  than  a  body  length,  a  front  is  a  reasonable 

(12) 

approximation.  Detonation  is  a  heat  front  of  a  very  special  nature. 
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A  heat  front  with  a  subsonic  leading  edge  located  next  to  a  wall  will  be 
preceded  by  an  oblique  shock.  The  oblique  shock  turn!*  the  flow  away  from  the  wall, 
anil  the  heat  front  turns  it  back  parallel  to  the  wall.  Recall  that  for  a  subsonic 
leading  edge  the  heat  front  can  influence  the  flow  upstream.  The  corresponding  changes 
in  pressure  are  shown.  The  net  change  in  pressure  is  positive  although  the  heat 
front  decreases  pressure. 

A  heat  front  with  a  supersonic  leading  edge  will  be  followed  by  an  expansion 
fan  if  the  front  is  adjacent  to  a  straight  wall.  Once  again  the  net  pressure 
change  is  positive.  The  heat  front  turns  the  flow  away  from  the  wall,  and  the 
fan  redirects  it  so  as  to  be  parallel  to  the  wall. 

Adding  heat  to  a  supersonic  stream  decreases  the  Mach  number.  When  sufficient 
heat  has  been  added,  the  normal  component  of  the  Mach  number  downstream  of  the 
heat  front,  M^,  will  become  unity.  The  amount  of  heat  required  to  make  s  1  is 
known  as  "critical"  amount  of  heat.  With  a  heat  front  having  a  supersonic  leading 
edge,  heat  can  be  added  until  =1.  If  more  heat  is  added,  the  flow  illustrated  in 
Figurel6(b)  changes  to  that  shown  in  Figure  16(a).  The  oblique  shock  decreases 
the  flow  Mach  number  changing  the  leading  edge  of  the  heat  front  from  supersonic  to 
subsonic. 

Adding  heat  to  a  subsonic  flow  or  a  heat  front  with  a  subsonic  leading  edge 
modifies  the  flow  upstream.  One  cannot  specify  the  flow  immediately  ahead  of  the 
front.  Critical  heat  addition  does  not  have  the  same  impact  in  the  subsonic  case 
since  the  flow  upstream  will  be  modified  to  adjust  to  the  large  beat  addition. 

•  A  detonation  wave  ear.  bo  thought  of  as  a  shock  wave  followed  by  a  combustion  zone 
or  heat  front.  In  Figure.!#a)  if  the  oblique  shock  and  the  heat  front  had  the  dame 
angle,  it  would  be  an  oblique  detonation  wave  provided  the  Mach  number  and  amount  of 
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heat  reieaued  were  correct. 
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FIGURE  16.  HEAT  FRONTS  AND  COMBUSTION  FANS. 
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In  the  Prandtl -Meyer  expansion,  conditions  axe  constant  along  a 
characteristic.  Centered  combustion  fans  are  possible  with  conditions 
constant  along  a  ray.  However,  the  ray  i3  not  a  Mach  line,  and  the 
normal  component  of  Mach  number  need  not  be  unity.  Combustion  fans  could 
occur  in  subsonic  as  well  as  supersonic  flow. 


If  Mr,  the  normal  component  of  Mach  number,  exceeds  unity,  the  fan  will 
turn  the  flow  away  from  the  wall  like  a  compression.  This  is  illustrated  in 
Figure  16(c)  on  the  left-hand  side.  If  Mn  is  less  than  unity,  the  flow  would 
be  turned  downward  as  shown  in  Figure  16(c)  on  the  right-hand  side. 

Starting  with  the  continuity  equation  in  cylindrical  coordinates  (r,§,z) 
and  setting  8(  )/8z  =  0  and  8(  )/5r  -  0,  one  can  derive  the  continuity 
equation  for  a  combustion  fan.  It  is 


The  velocity  components  are  illustrated  in  Figure  17.  Heat  addition  does  not 


cause  an  irrotational  flow  to  become  rotational.  One  additional  equation  is 


the  condition  for  irrotational  flow  which  is 


(57) 


The  rate  of  heat  release  controls  d^np/df.  One  simple  combustion  fan  occurs 
for  d-tnp/dp  =  constant.  For  this  case,  Equations  (56)  and  (57)  can  be 


readily  integrated. 


There  is  a  particularly  interesting  expression  that  can  be  derived  for 
the  pressure  coefficient  associated  with  a  heat  front  in  hypersonic  flow. 
Consider  a  heat  front  at  angle  9  relative  to  the  main  stream.  The  velocity 
components  and  angles  are  illustrated  in  Figure  IS.  The  momentum  equation  is 
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2  2 

?!  +  P^  =  P2  +  P2u2  (58) 

which  can  be  rearranged,  since  from  continuity 

p2  “  “  plUl  -  P2U2  = 

From  the  geometry  illustrated  in  Figure  18 

V2  _  ”l  ~  u2 
sinj^  -  ©j  sin  6 

Noting  that  sin  0  and  combining  Equations  (59)  and  (60)  leads  to 

P2  -  P-^  =  P^w2  'tan  6  sinS  (61) 

If  =  M2  sin2©  »  1,  heat  addition  has  little  influence  on  w  or  p.  Hence 
W1  -  w2*  ^or  hyP0rsonic  Equation  (61)  becomes 

p2  “  P1  =  2(lPlWl)tan  6  sin& 

and  the  pressure  coefficient  is 

0^  =  2  tan  9  sinfi  (62) 

This  can  be  compared  with  the  Newtonian  expression  which  is 

c  =  2  sin2 6  (6?) 

P 

Equations  (62)  and  (63)  provide  sufficient  tools  to  explore  the  external 
burning  ramjet  shown  in  Figure  19.  Angle  5  is  shown  in  the  top  figure  and 
applies  to  Equation  (63). 

Heat  Addition  and  Vorticity 

In  this  section  it  is  shown  that  heat  addition  does  not  change  an 
initially  irrotational  flow  to  a  rotational  flow.  The  curl  of  the  velocity 
vector  should  be  zero.  For  subsonic  flow,  one  uses  the  equations  from  Tsien 


*1*1  *  p2U2 

P1U1*U1  -  u2)  (59) 

(60) 
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and  Beilock^^ ;  and  for  supersonic  flow,  one  uses  Equations  (22)  and  (23). 

curl  7=  e^(f~'  "  fy~) 

Substituting  Equations  22  and  23  into  Equation  (64)  yields 


curl  v  =  e 


i|j*2yp^[~  2n^x  ”  $y)  (x  -  py)]  +  ^2ypg^  2n^x  "  ^  (*  “  fyjj} 


Recall  that  p2  =  M2  -  1.  In  the  above  analysis,  5  ha3  been  represented  by 


&(x  -  py)  =  axp[-  n(x  -  py)2]  (66) 

Since  curl  v  vanishes,  one  concludes  that  heat  addition  does  not  generate 
vorticity. 

Mass  Balance  with  Heat  Addition 

To  gain  additional  insight  to  diabatic  flows,  it  is  worthwhile  to  examine 
some  simple  flows  using  familiar  methods.  One  such  case  i3  a  mass  balance. 

Consider  a  heat  addition  region  as  shown  in  Figure  20  with  a  control  volume. 

Half  of  the  control  volume  can  be  used  for  the  mass  balance.  This  half  has 
width  w.  Writing  a  continuity  equation  for  mass  flux  across  control  volume 
yields 

P*u*w  s  (p  +  p')v’(L  -  ft)  +  2(p  +  p'Jv'p-t  +  po#ue0(w  -  1/2)  +  ( p o0  +  p^)u  1/2  (67) 

where  p'  and  v'  are  averages  of  perturbation  quantities  along  p l  illustrated 
in  Figure  20.  The  density  perturbation  in  the  wake  is  p*.  Equation  (67) 
reduces  to 

pv'L  =  -  phiJ/2  (68) 

when  second  order  terms  (e.g. ,  p'v')  are  neglected.  This  can  be  rewritten  as 
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A 


p  21 


Substitution  of  Equations  (40)  and  (42)  in  (69)  gives  an  equality,  as  one 


expects. 


EXTERNAL  BURNING;  THRUST  MINUS  DRAG 


Drag — Subsonic  and  Supersonic 

In  subsonic  flow  of  a  perfect  fluid  around  a  body,  there  is  no  force  on 
the  body,  either  lift  or  drag.  Motion  of  a  body  through  a  fluid  would  require 
work  if  there  were  a  drag.  This  energy  would  appear  in  the  fluid.  For  a 
perfect  fluid,  the  perturbations  in  the  fluid  decay  rapidly  at  large 
distances  from  the  body.  If  the  body  has  been  moving  for  a  long  time  and  if 
there  were  drag,  there  should  be  significant  motion  of  the  fluid  at  large 
distances. 

In  supersonic  flow,  there  is  significant  motion  at  large  distances  even  in 
a  perfect  fluid.  The  waves  caused  by  body  motion  move  out  from  the  body  eventually 
decaying  to  sound  waves.  Far  from  the  body,  the  flow  consists  of  outgoing  sound 
waves.  These  waves  carry  away  energy.  Work  overcoming  drag  su  j  ilies  the 
energy. 

Relation  Between  Flux  of  Entropy  and  Drag 

It  is  well  known  (See,  for  example,  K.  Oswatitsch.  )  that  the  flux  of 
entropy  through  a  control  surface  can  be  related  to  drag  of  a  body  within  the 
control  surface.  To  get  a  feeling  for  the  extent  of  the  waves  causing  drag,  a 
sample  problem  was  worked.  The  results  are  shown  in  Figure  21. 

Consider  a  20°  wedge  at  Mach  2,  as  shown  in  Figure  21.  The  wave  geometry 
is  shown  and  has  been  calculated  using  the  method  of  finite  waves.  Width  of 
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the  base  of  the  wedge  is  b.  The  lateral  distance  to  a  point  on  the  3hock 
wave  is  t.  Distance  i3  illustrated.  The  flux  of  entropy  in  and  out  of 
the  control  volume  (neglecting  the  wake)  is  equal  to  entropy  jump  across 
the  shock  wave  times  the  mass  flux.  The  first  expansion  wave  intersects  the 
shock  wave  at  point  2.  The  \n  re  between  points  1  and  2  accounts  for  25.2 
per  cent  of  wedge  drag.  The  shock  wave  between  points  1  and  1  accounts  for 
39.3  per  cent  of  the  wedge  pressure  drag.  It  is  necessary  to  go  out  179 
base  widths  to  obtain  99  per  cent  of  the  drag. 

The  conclusion  from  this  little  study  is  that  the  waves  are  important  to 
great  distances  from  the  body.  In  axisymmetric  flow,  the  bow  wave  would  decay 
more  rapidly. 

Linearized  Drag  Due  to  Waves 

A  wave  is  a  mechanism  whereby  energy  can  be  propagated.  The  energy  takes 
two  forms  in  the  linearized  or  acoustic  case.  These  are  flow  work  and  kinetic 
energy.  It  can  be  shown^^  that  an  integral  over  space  of  the  flow  work  equals 
the  integral  over  space  of  kinetic  energy.  As  an  example,  consider  a  half¬ 
diamond  airfoil  shown  in  Figure  22.  The  airfoil  i3  moved  exactly  one  chord 
length  from  position  ABC  to  position  DEA,  and  the  undisturbed  fluid  is  at  rest. 
This  motion  does  work  Dc  where  D  is  the  drag  and  c  is  the  chord.  The  work  appears 
as  energy  within  the  waves.  Changes  in  the  fluid  motion  occur  normal  to  the 
waves.  The  disturbance  originating  at  A  now  appears  at  A1,  B  at  B’,  C  at  C,  etc. 
The  motion  within  the  wave  below  A’B’C’  is  due  to  earlier  airfoil  motion  to  the 
right  of  position  ABC.  The  work  done  by  drag  moving  one  chord  length  appears 
within  the  volume,  V,  bounded  by  A'B'C'DKA.  Consequently 
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where  u  and  v  are  perturbation  velocities.  When  there  is  heat  addition,  as 
illustrated  in  the  lower  part  of  Figure  22,  Equation  ('’O)  changes.  The  force 
on  the  body  is  now  thrust  minus  drag.  Heat  has  been  added.  There  is  a 
left-running  wave,  AG,  which  moves  into  the  fluid  above  the  airfoil.  Energy 
is  added  to  the  fluid  bounded  by  ACG.  Heat  is  added  in  a  volume  bounded  by  BFC. 
Downstream  of  this  region  there  is  a  wake.  Note  that  in  the  region  EAC’B1  the 
mass  sinks  cause  u1  and  v*  to  be  positive*  whereas  the  heat  sources  cause  u" 
and  v"  to  be  negative.  In  region  EAC'B'  the  mass  sinks  decrease  p  whereas  heat 
increases  p.  The  energy  equation  associated  with  the  waves  now  is 


r 

(T  -  D)c  +  (heat  added)  =  ^  pl(U  +  u*  ♦  u")2  +  (v*  +  v")2JdV 

J 


V 


pd(~)  j  ( energy  )dV  + 
'WAKE 


(7l) 


( energy  )dV 


AGO 


Equation  (71)  is  more  conceptual  than  quantitative.  There  are  some  questions 
yet  to  be  answered.  Does  the  integral  of  flow  work  equal  the  integral  of 
kinetic  energy  for  waves  caused  by  heat  addition?  Incorporating  heat  addition 
into  the  waves  increases  complexity.  Can  this  method  of  drag  and  thrust 
-  x  ation  lead  to  useful  results  and  new  insight?  These  questions  will  be 


*We  U3e  the  notation  that  masssources  cause  perturbation  velocities  u’  and  v1; 
heat  sources  cause  perturbation  velocities  u"  and  v'1. 


considered  in  the  section  Wave  Energy  Due  to  Heat  Addition  which  follows 
shortly. 

Linearized  Momentum  Relations 

In  two  dimensions,  x  and  z,  the  drag  formula  in  three  dimensions  is 
simplified.  It  becomes 

D  =  -  (* (p  -  P^)dz  -  pCn  ♦  u)£(U  +  u)dz  -  wdxj 

c  Jc 

The  corresponding  lift  equation  is 

l'  j(p '  p»)dIt  *  pwjwdx  -  (U  ♦  u)dzj 
C  C 

where  C  denotes  a  closed  contour  surrounding  the  vehicle  or  wing.  For 
small  perturbations  due  to  mass  sources  and  sinks,  the  above  equations 
become 


D-  P, 


,  1  ro2  2  ^  21 

uwdx  +  u  +  w 

*  n  n 


2  2 

Note  that  drag  is  due  to  second  order  terms,  i.e.,  uw,  u  ,  and  w  .  First 
order  terms  drop  out.  For  lift  the  small  perturbation  form  is 


{udx  ♦  vdz)  =  paUT 


In  a  following  section  the  equivalent  formulas  for  L  and  D  will  be  derived 
for  combined  mass  and  heat  addition.  First,  consider  drag  or  thrust  on  a  heat 
addition  zone. 


To  illustrate  the  application  of  the  heat  addition  equations,  let  us 
calculate  the  thrust  of  a  line  source  of  length,  t,  oriented  parallel  and 
normal  to  the  stream.  Equation  {^6)  is  taken  from  Chapter  7  of  Liepmann 
and  Roshko;^^ 


D  =  - 


* 

(PuX)nkukdA  “ 

JA 


J  pr^dA 

A 


(76) 


In  their  notation  =  x  and  x0  =  y.  The  control  volume  is  shown  in  Figure  23. 
If  Equation  (^6)  is  negative,  one  gets  a  thrust.  Using  the  numbers  for 
surfaces  in  control  volume: 


D  =  - 


(p  U)(-  U)dA 


-I 


(e«>- 


pJME)U2dA 


|  (0  -  ■•)(,.  -  1)(-  V)dA  -  (77) 

*5 

- 

<U  -  u*)(Poo  "  P'wAVE^+  1)^v’)dA 

A6 

The  integral  over  A.^  can  be  integrated.  The  integral  over  A?  can  be  split 
into  two  parts,  and  the  part  with  pw  can  be  integrated.  Equations  (22)  through 
(25)  were  derived  dropping  terms  in  second  order  in  v'/U,  etc.  The  integrands 
of  A*,  and  A^  reduce  to  p^Uv'.  Consequently,  Equation (71/)  becomes s 
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(78) 


We  now  need  expressions  for  P'y^g  and  v’»  which  we  obtain  from 
Equations  (25)  and  (23)  respectively.  Proceeding  with  p' 

l 


PWAKE^X,y^  = 


.(l*.  .7-  ill  5(y  _  yt  )l(x  -  x’Jds 

a  U 


(79) 


0 


where  s  is  distance  along  the  source;  x  and  y,  the  coordinates  where  you  want 
to  evaluate  Py^ggi  and  x’  and  y',  the  coordinates  of  an  element,  ds,  of  heat 
source.  When  x>  x’,  1*1;  since  we  are  evaluating  A, -integral,  x  will  be 
greater  than  x'.  Also  ds  =  dx‘  for  the  case  at  hand. 

,1 


iiinr/y)  =  6(y  -  y') 
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dx»  =  g(y  .  y.) 

a  U 


(80) 


Substituting  Equation  (80)  into  A,-integral  of  Equation  (78),  we  get 


p,WAKEU  dA  s'' 


A, 


2.. 

a  U 


(81) 


3  ? 

In  a  similar  manner,  an  equation  can  be  obtained  for  v*  from  Equation  (23): 


f' 


v'(x,y)  *  (  6[(x  -  x«)  -  E(y  -  y«)]ds  (82) 

yj 

0 


Substituting  Equation  (82)  in  A,.-integral  of  Equation  (^8)  yields 


57 


r  rl 

6[(x-x«)  -  p(y  -  y*)jdx*dx 

J  J 


a5  0 

The  coordinate  system  is  chosen  so  that  the  heat  source  starts  at  x'  *  0  and 
ends  at  x*  =  -t.  Also  y'  =  0.  The  coordinate  of  point  c  in  Figure  21  is 
(py^., y^)»  and  the  coordinate  of  point  d  is  'py,.  +  -t,yj).  Consequently,  the 
limits  of  integration  are  as  shown 


D<!^2  J  J6[(x-X')-  fy5] 


dx’dx 


A  change  of  variable  helps  evaluate  the  integral;  let  z  =  x  -  py^. 

I  l 

pyxwh  r  «<»  -  *,>dzdx'  - 

%  % 


0  0 


The  integral  in  Equation  (19)  is  in  the  standard  form  as  shown  in  the 
Appendix  so  that  it  is  readily  apparent  that  it  has  value  l.  Now  combine 
Equations  (78),  (81),  ard  (85)  to  give 


D  =  ^  ~  2  (86) 

a  \  2YP  / 

which  vanishes  since  a  =  Yp/p*  Supersonic  heat  addition  in  an  unconstrained 
(no  walls)  flow  gives  zero  thrust  or  drag. 

Now  look  at  the  case  where  the  heat  source  is  of  length  l  but  normal  to 
the  streamlines.  This  is  illustrated  in  Figure  24.  Equations  (78),  (79),  and 
(82)  apply;  however,  ds  is  different  as  is  the  limit  of  integration.  Put 
the  origin  of  coordinates  at  the  center  of  the  heat  source.  Equation  (79)  becomes 
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(87) 
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1/2 

r^3  I  6  (y  -  y')dy’ 

u  J 


-1/2 


The  A^-integral  in  Equation  (78)  becomes 

,4/2 


I"2 


'i2  -!X-j-iia  6(y  -  y')dy’dy  =  <*  -  ^1to 

a4U  J  a 

'A,  -  4/2 

since  the  integrations  over  y'  and  y  overlap  over  the  region  -  1/2  to  1/2. 
For  the  perturbation  velocity,  the  equation  for  v*  becomes 


(88) 


.V2 


v'(x,y)  = 


-  4/2 


"  ^y5  **  y,5dy1 


(89) 


where  y^  is  the  coordinate  for  plane  A,..  The  A^-integral  is 

‘4/2 


PooUv«dy  =  I  /  ~^qP<"U6[x  -  p(y$  -  y*)]dy'dx 


(90) 


The  limits  of  integration  for  x  are  g(y^  -  4/2)  to  {i(y*.  «■  4/2).  To  obtain  a 
standard  form  for  the  integrand,  we  must  eliminate  £  in  the  6 -function.  This 
we  do  by  letting  z  *  fS(y*j  -  y' )  and  dz  =  -  fdy*.  Equation  (90)  is  now 
r  r 


^  2YPP 


5  (x  -  z)dzdx  =  p(y5  +  §)  -  p(y5  -  \)  (91) 
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The  limits  of  integration  for  both  a  and  x  in  Equation  (91)  are  f(y^  -  1/2) 
to  P(yc  +  V2).  Combining  Equations  (78),  (88),  and  (91),  one  finds 


(Y-l)qlu  /(y  ^qtP<8U\ 
D  *  - 5 - 2 


0 


(92) 


a  \  2yp 

which  is  what  we  expected. 

It  is  interesting  to  draw  the  streamlines  in  the  vicinity  of  the  heat 
sources;  see  Figure  25.  For  the  parallel  heat  source  there  is  a  conceptual 
difficulty  which  has  a  resolution.  The  "stagnation"  streamline*  C  goes 
straight  back  since  v*  is  an  odd  function  and  must  be  zero  at  y  =  0.  One 
can  think  of  a  streamtube  bounded  by  streamlines  B  and  D.  Continuing,  one 
can  shrink  this  streamtube  so  that  there  is  very  little  mass  flow  through  it. 

Yet  the  streamtube  diverges  at  angle  9.  A  finite  amount  of  heat  is  added  to  an 
infinitesimal  stream  yielding  a  finite  angle  8.  Strange.*  One  avoids  this 
difficulty  by  using  volume  heat  addition.  The  case  of  the  heat  source  normal 
to  the  flow  is  more  easily  digested.  Streamlines  F  and  J  are  acted  upon  by 
characteristics  of  a  single  family.  Streamline  G  passes  through  the  heat 
source  into  a  region  where  it  is  influenced  by  both  left-running  (gives  +  v*) 
and  right-running  (gives  -  v*)  characteristics.  In  this  region  the  streamline 
is  not  deviated  from  its  original  direction.  Streamline  G  then  moves  into  a 
region  where  only  left-running  characteristics  act  on  it.  It  is  deflected  by  an 
amount  0.  Finally  G  moves  into  the  wake  where  it  accelerates  and  assumes  its 
original  direction.  The  speed  along  G  is  (0,  U  -  2u*,  U  -  u*,  U)  33  it  moves 
from  region  to  region.  Streamline  H  is  in  the  plane  of  symmetry;  as  a  result, 
it  does  not  deviate  from  freestream  direction. 


*A  body  will  stagnate  the  flow;  however  a  heat  source  does  not  do  so. 
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Wave  Energy  Due  to  Heat  Addition 

The  model  for  determining  wave  energy  due  to  heat  addition  is  shown 
in  Figure  26.  A  heat  addition  zone  EFGH  is  translated  its  own  Length  L  in 
time  L/ua.  Its  new  position  is  AFEL.  The  heat  released  is  (btLHL/u^). 

Since  there  is  no  force  on  the  heat  addition  region,  no  work  is  done  moving  it 
from  EFGH  to  AFEL,  The  question  of  forces  on  heat  addition  zones  was 
discussed  in  the  preceding  section.  All  of  the  heat  energy  released  must 
appear  in  the  waves  and  in  the  wake. 

In  the  time  L/u,,,,  the  disturbance  at  J  has  moved  to  point  K.  The 
disturbance  at  F  has  moved  to  B.  The  dashed  lines  AG  and  BD  would  be  stream¬ 
lines  if  EFGH  were  stationary  and  tne  fluid  moved  by.  The  energy  content  of 
space  bounded  by  ANCEL  is  the  same  as  that  inside  FJKHE.  Part  of  the  energy 
added  appears  in  the  region  ABDC.  This  is  the  wave  energy.  The  remainder 
appears  in  the  wake  in  the  region  EFGH.  From  B  to  P,  work  is  done  on  the 
fluid  to  compress  it.  From  K  to  D,  work  is  done  by  the  fluid  as  it  expands. 
The  flow  work  cancels.  We  are  now  ready  to  write  the  energy  equation. 


hlL2  s  2  (  kinetic  energy\  /  energy  in 
u  in  ABDC  wake  EFGH 


(93) 


The  factor  2  comes  from  the  fact  that  waves  go  both  up  and  down  in  Figure  26. 
The  kinetic  energy  in  ABDC  is 
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(94) 


where  dx  and  dy  have  been  replaced  by  a  coordinate  system  x  =  £  and  y  =  sin 
The  symbol  u^  represents  the  perturbation  velocity  normal  to  the  wave.  It  is 


eaual  to 


63 


U2  =  U*2  +  v'2  (95) 

n 

where  u'  and  v*  are  given  by  Equations  (39)  and  (40)  respectively.  The 
magnitude  of  v*  has  been  plotted  in  the  upper  portion  of  Figure  20.  Inserting 
Equation  (95)  into  Equation  (94)  and  recognizing  the  functional  form  of  v' 

(and  u* )  as  shown  in  Figure  20,  one  arrives  at 


2KE  in  ABDC 
Lp^sin  y.  cos  ji 


Pf.  2  ,  2,  (■ 

)P)  t(^)Jj2dr*(U»2*V»2)  j  W 


The  value  of  v*  has  been  indicated  in  Figure  20.  It  is  the  maximum 
m 

perturbation  velocity.  Integrating  Equation  (96)  yields 


2KS  in  ABDC 
tLp  sin  y  cos  }i 


-)  -  f(ui2  *  *  k2  ♦  7i2)(L  -  & 


The  energy  in  the  wake  i3 


(wake^Gh)  "  po*CpL^Tw  "  T*^  "  p<wcpT®L^Te(> 


The  distance  FG  is  the  length  L.  According  to  small  perturbation  results, 

the  wake  has  freestream  velocity  and  pressure  but  reduced  density  and  increased 

temperature.  Knowing  these  facts  and  using  p  =  pRT,  Equation  (98)  becomes 

/energy  in\  -  Lj  ,  (go) 

^wake  EFGHj  pTaLipw 

The  density  perturbation  in  the  j^ake  p^  can  be  obtained  from  Equation  (44). 

Substituting  Equations  (39)  and  (40)  into  Equation  (97),  substituting 
Equation  (44)  into  Equation  (99),  and  combining  the  results  in  Equation  (93) 
yields 
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WAKE  ENERGY 


From  Equation  (100 )  it  is  apparent  that  the  energy  appearing  in  the  waves  is 


2  2 

second  order.  The  wave  energy  was  based  on  u'  +  v*  which  is  second  order. 

mm 


The  wake  energy  reduces  to  QlL  /u^.  If  Equation  (100)  is  to  be  an  equation, 
the  wave  energy  must  be  neglected. 

There  are  some  important  conclusions  from  this  exercise.  First,  the 
wave  energy  i3  second  order.  Second,  most  of  the  energy  due  to  heat  addition 
appears  in  the  wake.  Ehergy  radiated  away  by  waves  is  small  compared  to 
energy  deposited  in  the  wake. 

If  one  examines  the  linearized  momentum  equations  in  two  dimensions,  as 
given  in  the  section,  Linearized  Momentum  Equations,  one  finds  drag  is  related 
to  an  integral  of  second  order  terms.  As  shown  in  the  discussion  of  Figure  21, 
drag  can  be  related  to  the  waves  generated  by  motion  of  a  body.  A  conclusion 
is  that  wave  energy  resulting  from  either  drag  or  heat  addition  is  second  order 
in  the  perturbation  quantities. 

Control  Volume  Approach  to  Thrust  Minus  Drag  with  Heat  Addition 

As  an  example  of  a  two-dimensional  planar  flow,  examine  Figure  27.  A 
half  diamond  has  a  region  of  volume  heat  addition  indicated  as  a  shaded  area. 

The  dotted  lines  outline  an  image  heat  addition  region  above  the  shaded  region. 
Part  of  a  control  volume  boundary  is  shown  as  a  dashed  line.  Mach  lines  are 
drawn  from  the  body  to  the  boundary  of  the  control  surface.  Below  the  boundary 
are  plotted  the  perturbation  quantities  for  that  location  on  the  boundary.  At 
the  very  bottom  are  the  pressure,  source  strength,  u*  and  v’.  These  are  constant 
to  mid  chord  and  then  switch  to  values  opposite  in  sign. 
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FIGURE  27.  CONTROL  VOLUME  APPROACH  TO  DRAG  AND  THRUST. 


The  perturbation  quantities  due  to  heat  addition  are  plotted  about  the  right- 
hand  side  of  mid-page.  Heat  addition  increases  pressure,  slows  the  flow,  and 
deflects  it  downward.  If  the  two  plots  were  to  scale,  the  positive  pressure 
increment  due  to  heat  addition  would  exceed  the  negative  increment  due  to  flow 
expansion  (mass  sinks).  The  body  would  have  positive  thrust  minus  drag. 

The  magnitudes  of  the  perturbations  due  to  heat  addition  are  proportional  to 
the  length  of  the  characteristic  (Mach  line)  passing  through  the  heat  addition 
regions  both  real  and  image.  Note  that  p',  p”,  u',  u”,  v’,  and  vrt  are  of 
opposite  signs. 

Let’s  look  at  drag  and  lift  contributions  due  to  the  portion  of  control 
volume  shown.  In  Equation  (72)  dz  is  zero.  Hence  the  drag  contribution 

D  *  J"  p(U  +  u’  +  u" )(w'  +  w”  )dx  (101 ) 

Multiplying  the  terms  yields 

D  =  J  (pUw*  ♦  pUw”  +  pu'w’  ♦  pu’w"  ♦  pu"w’  *  purtw")dx  (102) 

1  2  3  4  5  6 

The  following  summary  gives  the  contributions  (+),  0,  or  (-). 

Term  I  2  3  4  5  T~ 

0  to  c/2  (-)  0  (+)  0  0  0 

c/2  to  c  (+)  (*  '  (+)  (-)  (-)  (■»•) 

The  integral  of  pUw’  over  a  full  chord  gives  zero;  however,  the  integral  of 
pUw"  makes  a  large  thrust  contribution.  Terms  3  and  6  give  drag,  whereas  terms 
2,  4,  and  5  give  thrust. 


Lift  is  obtained  from  Equation  (7?),  Once  again  for  the  segment  of 
control  volume  illustrated,  dz  is  zero.  The  lift  contribution  becomes 

L  =  (p '  +•  p")dx  ■*■  I  p(w'2  +  2w'w"  +  w"2)dx  (103) 

J  1  2  3  4  5 

The  following  summary  gives  the  contributions  of  the  various  terms. 

Term  1  2  3  4  5 

0  to  c/2  (+)  0  (♦)  0  0 

c/2  to  c  (-)  (♦)  (»)  (-)  (♦) 

Tne  pressure  p'  averages  to  zero,  whereas  prt  gives  a  lift.  The  cross  term 
u'w”  reduces  lift. 

To  optimize  the  heat  addition,  the  airfoil  shape  must  be  expressed 
analytically.  H  ,at  addition  also  needs  to  be  described  by  some  function  of 
x  and  z.  Heat  rolease  need  not  necessarily  be  uniform. 

Thrust  from  Cycle  Point-of-View 

The  thrust  due  to  a  propulsion  device  is  given  by 

F  =  m(w  -  vm)  *■  AQ(pe  -  p^)  (104) 

If  one  evaluates  the  equation  at  the  Tret f tzplane,  the  pressure  term  may  or 
may  not  drop  out.  This  is  ono  item  which  Oswatitsch^  ^  discusses.  In 
’Equation  (104),  F  is  thrust, and  m  is  raaso  flow  rate  of  fluid  which  receives 
a  velocity  increase  w  -  over  the  freostream  velocity  w^ 

Consider  three  flows  as  shcun  in  Figure  28.  Heat  is  added  at  stagnation 
conditions  and  at  some  Mach  number  other  t:mn  zero.  The  pressures  at  the 
exhausts  are  all  equal,  i. e.,  pro=  p^  -  =  p!>.  The  Mollier  Chart  which  depicts 
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ISENTROPIC  WITH  NO  HEAT 
ADOmON;THIS  IS  A  REFERENCE 


(b)  ISENTROPIC  WITH  HEAT 
ADDITION  AT  STAGNATION 


(c)  ISENTROPIC  TO  STATION  2;  ADD 
HEAT  2  T03;  ISENTROPIC  3  T04. 


FIGURE  28.  THREE  FLOWS  WITH  TWO  DIFFERENT  WAYS  OF  ADDING  HEAT. 


the  three  flows  is  shown  in  Figure  29.  It  is  apparent  from  energy  considerations 


Continuing 


dq  *  ¥2  =  2W4  *  (h4  “  h2} 


(105) 


2^w4  *  v2~^  =  wdw  =  dq  -  (h  -  hp 


(106) 


dw  >  0 


V  if  dq  >  (h.  -  h.*)  =  dh 
w,  >  w'  4 


The  slope  of  a  p  =  constant  curve  is 


Oh]  _ 
Os 

n 


(1a7) 


(103) 


Combining  Equations  (106)  and  (108)  leads  to 


,  ap.  \ 

wuw  =  dq  -  dh  =  dq  -  ■j^-l  ds 


(109) 


The  change  in  entropy  of  the  flow  is  given  by 


ds  =  fa 

q 


(110) 


A  subscript  q  is  used  to  indicate  this  i3  the  temperature  at  which  neat  is 
added.  Equations  (109)  and  (110)  yield 


wdw  =  dq(l  -  ~) 


(111) 


If  one  compares  flow  (a)  with  flow  (b)  of  Figure  28,  then  the  temperature 
ratio  in  Equation  (113)  is  T*/T,p.  If  one  compares  How  (a)  with  flow  (c),  the 
temperature  ratio  is  T'/'d.  To  have  dw  >  0,  one  must  have  T  <  T  . 

<.  n 
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As  a  numerical  example,  consider  flight  at  M  *  0.9.  Heat  is  added  at 
M  =  0  in  flow  (b)  and  at  M,,  =  0. 5  in  flow  (c).  Enough  heat  is  added  to 
accelerate  the  flow  so  ■‘•hat  =  0.6.  The  same  amount  of  heat  i3  added  in 
flow  (b).  The  results  that  w^/w ^  =  1.055  and  wlj/w],  =  1.091. 

Equation  (ill)  is  not  usually  applicable  to  supersonic  or  hypersonic 
flow.  The  entropy  change  in  the  exhaust  was  attributed  solely  to  heat 
addition.  In  supersonic  flow,  shock  waves  cause  entropy  changes.  Further¬ 
more  heat  addition  may  alter  the  shock  wave  geometry. 

Thrust  and  Drag  in  Various  Flow  Regimes 

For  subsonic  flow,  the  pressure  term  in  Equation  (104)  is  zero.  For 
supersonic  flow,  the  pressure  term  if  evaluated  in  the  Trefftzplane  is  not 
zero;  however,  it  is  extremely  small.  To  see  that  the  pressure  is  not  zero, 
examine  the  flow  illustrated  in  Figure  30.  The  entropy  in  the  Trefftzplane 
is  greater  than  ambient;  ds  >  0.  If  dp  were  zero,  then  pwdw  equals  -  Tds 
so  that  dw  <  0.  Entropy  is  related  to  p  and  T  by  Tds  =  c^dT  -  dp/ p  -  c^dT. 
Consequently  dT  >  0.  From  p  -  pRT,  it  follows  dp/p  <  0.  The  continuity 
equation  pwA  must  be  satisfied.  As  shown  in  Figure  30,  the  streamtube  area  A 
does  not  change  significantly.  The  changes  in  both  p  and  w  cannot  be  negative 
and  still  satisfy  continuity.  If  dp  >  0,  then  dw  continues  to  be  less  than 
zero  and  dT  >  0.  Also  dp/p  >  dT/T.  Consequently  dp  >  0. 

To  account  for  the  drag,  the  area  A  must  be  quite  large.  (See  Figure  21.  ) 
The  streamtube  influenced  by  a  propulsion  system  will  be  small  as  shown  in 
Figure  30.  It  is  an  area  of  approximately  Ao  in  Equation  (104).  The  pressure 
at  the  Trefftzplane  is  due  mainly  to  drag  and  is  spread  over  area  A,  albeit  not 
uniformly.  The  pressure  term  in  the  thrust  equation  is  applicable  to  area  Ag. 
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Since  thrust  equals -drag,  the  pressure  term  will  be  of  the  order  of  Ag/A.  It 

can  be  neglected  for  subsonic,  transonic,  and  supersonic  flow. 

For  slender  bodies  in  hypersonic  flow,  dw/w  is  small.  If  the  thickness 

ratio  is  T=  ( thickness )/( body  length),  then  dw/w  is  of  the  order  Heat 

addition  also  has  little  influence  on  dw/vr,  a  fact  that  can  be  verified  by 

examining  Figure  7.3  in  Volume  I  of  A.  H.  Shapiro's  The  Dynamics  and  Thermo - 

(7) 

dynamics  of  Compressible  Fluid  Flow.  In  hypersonic  flow,  heat  addition 

causes  large  changes  in  pressure  and  static  temperature. 

In  hypersonic  flow,  the  shock  waves  are  close  to  the  body.  The  drag  is 

concentrated  in  a  small  area,  A^;  see  Figure  30.  Thrust  pr dicing  devices 

would  influence  an  area  comparable  to  Aq.  Consequently,  the  pressure  term  in 

Equation  (104)  must  be  retained  as  a  significant  factor. 

One  can  calculate  forces  on  a  body  by  suitable  integration  of  pressure. 

Pressure  in  excess  of  ambient  on  forward-facing  surfaces  causes  drag;  define 

i  is  as  pressure  drag,  D  ,  Pressure  less  than  ambient  on  r<  "•'ward -facing 

P 

surfaces  causes  drag;  define  this  as  suction  drag,  D  .  The  ratio  D  /D  is 

s  s  p 

shown  in  Figure  31  as  a  function  of  Mach  number.  Suction  drag  is  very 

important  in  subsonic,  transonic,  and  supersonic  flow.  However,  as  M 

becomes  large  and  enters  the  hypersonic  region,  D  becomes  insignificant. 

s 

In  fact,  Newtonian  theory  assumes  ambient  pressure  p^  is  zero  with  the  result 

D  =  0. 
s 


A  propulsion  device  will  increase  pressure  on  rearward  surfaces  as  shown 
in  Figure  31.  Heat  addition  must  increase  pressure  so  that  the  area  under  the 
curve  labelled  "with  heat"  equals  the  area  under  the  curve,  D^.  In  supersonic 
flow,  the  thrust  device  simultaneously  cancels  a  large  Dg  and  balances  D^. 
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In  hypersonic,  flow,  the  thrust  device  does  th>->  same  thing  except  that  in 

hypersonic  flow  D  is  negligible.  Heat  addition  in  hypersonic  flow  balances 

D  but  counteracts  only  an  insignificant  D  .  This  has  implications  con- 
P  3 

cerning  propulsion  efficiency. 

One  more  point  is  that  in  hypersonic  flow  is  the  major  dr<ug  con¬ 
tribution,  and  this  can  be  related  to  the  bow  shock  wave.  In  supersonic  flow, 

D,  is  important  also  so  that  the  waves  originating  at  the  rear  end  of  the  body 
or  at  the  rearward-facing  s  'riiices  are  of  equal  importance  to  the  bow  shock 
wave  in  calculating  drag. 

Contrast  in  Thrust  at  Super jonic  and  Hypersonic  Speeds 

Consider  a  body  immersed  in  a  flow  as  shown  in  Figure  30.  The  force  on 
the  body  in  the  direction  of  the  freestream  is 

F  =  j  [(pu2  -  Pau^)  *  (p  -  p^)]dydz  (112) 

A 

If  there  is  no  energy  addition  and  if  A  is  the  area  of  the  large  streamtube, 

F  equals  D.  (it  has  been  assumed  that  expansion  waves  have  made  pressure 
increments  negligible  along  the  streamline  from  a  to  b  in  Figure  30. )  If  there 
is  energy  addition  to  the  stream  and  if  A  is,  once  again,  the  Trefftzpline  in 
Figure  30,  Equation  (ll2)  yields  thrust  minus  drag.  If  one  evaluates 
Equation  (112)  using  the  small  streamtube  in  Figure  30  with  flow  m,  then  F  is 
approximately  thrust.* 

*It  is  not,  of  course,  possible  to  separate  out  the  thrust  accurately  using 
Equation  (112).  If  it  were  possible,  how  much  more  simple  the  propulsion 
integration  problem  would  be. 
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The  continuity  equation  is,  neglecting  fuel  added, 


0  =  \  (pu  -  PeoU«)dydz 


(111) 


The  energy  equation  in  integral  form  is 


|pu(h  + 


2  2 
u  *■  v 


)  -  pwu«.(hc»  *  t  *  9)1  dydz 


(114) 


In  Equation  (114),  q  is  the  heat  added  per  unit  mass,  and  v  is  the  component 
of  velocity  normal  to  the  freestream.  By  combining  the  preceding  equations 
in  the  following  way 

u* 

(uM  x  thrust)  +  (hM  +  q  -  -^-)  x  continuity  -  energy  =  u^F 
one  obtains 


V  =  Pu^hoo  -  h  +  q)  ♦  u*(p  -  P.)  -  pu 


,  x2  a  2 

(u  -  uj  ♦  V 


- - 1  dydz 

2  J 


(115) 


Substituting  h  =  e  ♦  p/p  into  Equation  (115)  yields 


u/  -  <pu 


i[q  -  (e  -  ejJ  -  p(u  -  uj  -  p' 


(  2 
(u  -  u_;  +  v 


dydz  (116) 


There  are  three  terms  in  the  integrand  of  Equation  (116).  The  term  pu(e  -  e^) 
is  the  flux  of  energy  thrown  away  in  the  jet.  The  term  pu  [q  -  (e  -  ej]  is  the 
energy  converted  to  work  or  jet  kinetic  energy.  The  term  p(u  -  u.,)  is  the  flow 


work,  and  the  last  term  in  the  integrand  is  the  flux  of  kinetic  energy  of  the 
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jet  relative  to  the  stationary  surroundings.  The  left-hand  side  is  the  rate 

at  which  the  thrust  does  work  in  overcoming  drag. 

For  transonic  and  supersonic  flow,  all  terms  are  important.  In  hypersonic 

flow,  (u  -  O/ujbj  v/u®  <&  1.  Changes  in  velocity  due  either  to  shock  waves  or 

heat  addition  are  small.  However  Ap/p^,  AT/Tfc ,  or  Ap/p^  are  of  the  order  of 

unity  or  larger.  Based  on  these  comments,  one  can  neglect  the  jet  kinetic 

energy  in  hypersonic  flow.  In  hypersonic  flow,  the  tern  p(u  -  u^,)  varies  as 
2 

l/M  a3  '  ill  now  be  demonstrated.  Rewriting  the  term  as 

p(u  -  O  =  p«u«,(1  +  1} 

2 

one  can  introduce  the  Mach  number  dependence.  In  hypersonic  flow,  Ap/p  ^  K 

CD 

where  K  is  the  hypersonic  similarity  parameter  TM.  The  velocity  perturbation 
2 

varies  as  T  as  discussed  previc  .sly.  Substitution  of  these  quantities  into 
the  flow  work  term  yields 

p(u  -  uj  -  -^rd  +  k2)(k2)  (ii7) 

Compared  to  the  heat  addition  term  in  Equation  (116),  the  term  p(u  -  u,,,)  is 
small  for  K  1  and  K  =  0(l).*  Equation  (116)  reduces  to 


(118) 


which  is  valid  for  M  »  1  and  K  =  0(1).  In  hypersonic  flow,  Equation  (118)  tells 
us  that  the  thrust  is  simply  due  to  the  heat  added  less  the  amount  of  internal 
energy  thrown  away  in  the  exhaust.  In  supersonic  flow,  it  is  more  complex;  see 
Equation  (116). 


•Requiring  K  =  0(1)  is  not  a  severe  restriction. 
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APPLICATIONS  OF  EXTERNAL  BURNING 


K. 


I 


& 


if 


'1 


In  the  preceding  sections  some  of  the  analytical  tools  for  understanding 
external  burning  have  been  discussed.  It  is  worthwhile  to  look  at  some 
possible  applications  of  external  burning.  Four  applications  will  be  examined 
briefly  in  the  following  section;  these  are  forces  on  a  planar  airfoil  due  to 
heat  addition,  transonic  boattail  and  base  pressure  alleviation,  base  pressure 
modification  in  both  the  planar  and  axisymmetric  cases,  and  spin  recovery  of 
aircraft. 

Before  discussing  the  applications,  a  mode]  for  pressure  rise  due  to  heat 
addition  near  a  flat  plat  will  be  discussed.  This  mode]  uses  some  of  the 
results  from  EXTERNAL  BURNING;  BASIC  EQUATIONS  AND  SOLUTIONS  and  from 
EXTERNAL  BURNING;  THRUST  MINUS  DRAG.  The  model  i3  a  refinement  of  a  similar 
mode]  by  Billig. 

Model  for  Heat  Addition  Adjacent  to  a  Flat  Plate 

Figure  32  illustrates  the  geometry.  There  are  three  zones.  Zone  I  is 
the  freestream  which  is  supersonic.  Zone  II  is  the  heat  addition  region 
separated  from  done  I  by  a  heat  front.  Zone  II  is  bounded  by  the  wail,  heat 
front,  and  streamline  c  Zone  III  is  downstream  of  the  oblio.ue  shock  com¬ 
prising  the  region  below  streamline  c.  The  heat  addition  turns  the  flow  by 
an  angle  9.  The  streamline  deflection  angles  and  static  pressures  in  regions 

II  and  III  must  be  equal  along  streamline  c.  In  the  following  development, 

2 

terms  in  9  are  dropped  as  being  negligible. 

The  deflection  angle  is  approximately: 


x 


(119) 


The  continuity  equation  is 


p  V  -  p^y^  =  +  ~  ^2U2^2 

Across  the  oblique  shock  there  is  a  pressure  increase  equal  to 

P2  ~  Pi  Y^e 


(120) 


(121) 


Equation  (121)  defines  Define  Q’  as  the  heat  addition  per  unit  mass  of 
air  flowing  in  the  streamtube  defined  by  streamlines  a  and  c.  The  energy 


equation  is 


Vn  +  ^3  ”  CpTT3 


(122) 


or  as 


cPTT2  +  ^3  =  cpTTl 

(15) 

Using  Equation  (7.14)  from  Liepmann  and  Roshko,  one  can  derive  the 
following  momentum  equation 

-  P2V2  *  Wl  •  p2y2  *  p3y3  -  p2a‘  '  0 
Combining  equations  (119),  (120),  and  (123)  yields 

p2U2y2^U3  "  *  ^p3  ~  ♦  6x)  =  0 


(122) 


2  2 


(123) 


(124) 


A  solution  exists  for  u^  =  u2  and  p^  *  p2*  Let’s  examine  the  consequences 
of  that  solution.  From  the  continuity  equation  and  equation  of  state 

?2u2y2 


(125) 


Since  u^  =  u2  and  p^  =  p2>  this  reduces  to 


T3  _  y3 


(126) 
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The  energy  equation  become s 


12  12 
c  T~  ♦  tu~  +  Qi,  *  c  T,  ♦  ttu, 
p  2  2  2  21  p  3  2  3 


(127) 


c  T~  *  Qi,  »  c  T, 
p  2  '23  p  3 


since  u^  =  From  Equations  (126)  and  (127) 


y3  "  y2  Q23 


(128) 


In  the  notation  of  the  influence  coefficients  of  Table  IV,  Equation  (128) 


can  be  rewritten  as 


A  c  T 
P 


(129) 


In  view  of  the  equality  expressed  by  Equation  (129 ),  the  influence 
coefficients  give  the  result  dp/p  =  0  and  du/u  *  0. 


The  heat  added  by  the  heat  front  must  increase  pressure  in  the 


streamtube  a-c  by  an  amount  equal  to  pressure  rise  across  the  oblique  shock. 


Using  this  fact 


dp  yM^e  ytf*Q» 


*  Vl 


(130) 


Combining  Equations  (119)  and  (128),  the  heat  added  in  Zone  II  is 


Q23  *  CpT29x/jrl 


(m) 


The  difference  between  T^/T^  and  1  is  the  order  of  ©.  Hence 

«23  *  5pTlto/yl 


(132) 


■—  •  fa  n  .  ^ 


Total  heat  added  to  turn  the  flow  by  an  angle  0  is 

«5i  ■  «fc  *  «5i  -  9'1  -  f->  (!■>’) 

The  heat  addition  turns  the  flow  by  an  angle  0  and  increases  pressure 
by  the  amount  given  by  Equation  (130). 

Forces  on  &  Planar  Airfoil  in  Supersonic  Flight  Due  to  Heat  Addition 

A  problem  which  has  been  studied  by  several  investigators  is  external 

(17 ) 

burning  near  a  two-dimensional  wing.  One  such  study  is  that  of  Mager. 

Mager  uses  the  linearized  heat  addition  formulas,  Equations  (22)  to  (25),  to 
find  the  forces  on  an  airfoil  in  supersonic  flight  with  heat  addition  in 
adjacent  streamtubes. 

Consider  an  area  bounded  by  x.  and  x*  in  the  x-direction  and  by  a 
surface  and  a  parallel  side  at  distance  h^.  See  Figure  13.  Within  this  area, 
add  heat  Q.  From  Equation  (41 )  the  pressure  at  the  surface  is 


p' 


h  - 1  qm  \ 

2  a$p)  " 


sin  u  ds  * 


(y  -  l)QMh 


(114) 


P  \  2  a fPJ  J  r  ~  ”  a*p 

Where  does  the  factor  2  come  from?  There  must  be  an  image  heat  source  above  the 
surface;  otherwise,  there  would  be  a  v  with  flow  through  the  wall.  Combining 
Equation  (134)  and  the  definition  of  H  [h  is  defined  in  the  paragraph 
preceding  Equation  (11).]  gives 

p«  (y  -  l)HhM 
P  *  hi^xf  ”  xj_)a$P 

In  Equation  (135)  the  quantity  h  varies  from  0  at  x  *  to  h^  at  x  e  xi  +  ^i  co^  P* 
It  remains  constant  to  x  *  x^,  and  then  h  decreases  linearly  to  zero  at 
x  =  xf  +  \  cot  This  is  illustrated  in  Figure  34. 
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HEAT  H  UNIFORMLY 
RELEASED  IN  RECTANGLE 


FIGURE  33.  HEAT  ADDITION  ADJACENT  TO  A  SURFACE  IN  SUPERSONIC 


It  is  of  interest  to  examine  streamlines  near  the  surface  in  the  region 
of  heat  addition;  see  Figure  35.  The  streamline  passes  along  ABCD.  It  is 
shown  dashed  with  an  exaggerated  slope  below.  The  streamline  is  straight 
from  the  time  it  enters  the  heat  addition  region  until  it  arrives  at  A. 

This;  i3  so  since  aA  *  a' A.  From  A  to  B  the  right-running  characteristic  grows 
in  length  relative  to  aA.  Hence  the  streamline  is  curved  from  A  to  B.  From 
B  to  C  the  difference  in  length  b'B  -  bB  is  constant;  hence  the  streamline  is 


straight  but  at  an  angle  to  the  mainstream.  From  C  to  D  less  and  less  of  the 
right-running  characteristic  passes  through  the  image  heat  addition  region, 
and  v*  becomes  smaller  and  smaller  until  finally  the  streamline  is  horizontal 
at  D.  From  C  to  D  the  streamline  is  obviously  curved. 

When  the  streamline  i3  moved  closer  to  the  surface,  as  is  the  case  in 
the  lower  part  of  Figure  35,  the  deflection  is  considerably  less.  It  is 
obvious  thf streamline  on  the  surface  will  not  be  deflected  since 
aA  =  a'A  *  .  «  cC  «  c'C  *  d’B  *  dD  for  this  case. 

Gompa.  Tea  25  and  35;  'he  gains  some  insight  to  the  quandry 

discussed  in  >.  .lection  with  Figure  25.  When  the  heat  addition  is  confined  to 
a  mathematical  line,  as  in  the  top  of  Figure  25,  one  has  difficulty  because  of 
streamline  deflection  near  the  source.  However,  when  the  heat  is  distributed 
as  in  Figure  35,  there  is  no  problem. 

Based  on  Figure  34  and  Equation  (135),  it  is  apparent  that  one  has  in  hand 
the  tools  necessary  to  study  the  forces  on  bodies  due  to  heat  addition,  at 
least  in  the  linearized  case.  Refer  to  Mager's  paper  for  performance  of 


an  airfoil  with  heat  addition. 


Alleviation  of  Boattail  Drag  for  a  Plug  Nozzle 

In  the  Mach  number  range  0,8  ^  Mq  ^  1,2,  which  ii.  the  transonic  regime, 
the  nozzle  gross  thrust  coefficient  decreases  .r.if icantly.  A  factor  which 

may"  '*  this  transonic  dip  all  the  more  important  is  the  corresponding  increase 
in  tnt-  transonic  drag  coefficient.  Typical  data  for  a  plug  nozzle  are  shown  in 
Figure  ^6 ,  which  is  copied  from  the  NASA  Memorandum  by  Harrington. Part  (t) 
is  of  interest  here  while  pari  'b)  it  not  Two  flc;/  phenomena  cause  the  dip 
in  C^,  .  One  is  the  change  in  tn-ust.  figure  36(d)  showj  the  behavior 
of  plug  thrust.  Another  influx  is  4>uru*3  i-  rag  identified  in  Figure  36(c) 
as  primary  flap  drag.  Plug  thrusu  decree  and  boattail  drag  increases  in 
the  transonic  region. 

Is  there  anything  that  can  be  done  to  correct  the  loss  of  the  transonic 

V 

(19) 

Rabone  shows  that  by  going  to  small  plug  angles,  some  of  the  loss  of 

plug  thrust  can  be  recovered,  although  there  is  still  a  dip  near  M  =  1.2. 

(18) 

Harrington  uses  a  translating  shroud  that  has  a  pronounced  influence  on 
the  boattail  drag;  however,  the  loss  of  a  plug  thrust  partially  remains. 

An  alternate  approach  i3  to  burn  externally  on  the  boattail.  The  direct 
thrust  produced  by  high  pressure  on  the  boattail  hss  very  poor  SFC.  However, 
a  flow  interaction  may  occur  which  greatly  increases  thrust.  This  phenomenon 
ha3  been  termed  "wave  trapping”  by  Fuhs.  Figure  37  illustrates  some  essential 
features  of  the  interaction. 

With  heat  addition  on  the  boattail,  the  turning  angle,  9,  of  the  primary 
nozzle  flow  is  less.  See  the  angle  0  in  the  region  labelled  (l)  in  Figure  3*7. 

At  the  slipstream  between  the  external  flow  and  the  nozzle  flow,  the  expansion 

89 


waves  are  fleeted  as  expansion  waves  without  heat  addition.  These  reflected 
expansion  -  res  further  reduce  the  pressure  on  the  plug.  The  reflection  is 
indicated  point  (2)  in  Figure  ??,  With  heat  addition,  the  expansion  wave 
is  reflect  as  a  compression  wave  increasing  pressure  on  the  plug. 

In  re  >n  ("*)  without  heat  addition  the  pressure  is  low,  and  the  Mac:, 
number  is  -  -gs.  The  oblique  shock  (4)  has  a  shallow  angle.  It  is  necessary 
to  match  r>  .stream  pressure  and  flow  direction  for  both  the  external  and 
r.ozzle  flow  At  the  match  condition,  the  pressure  rise  across  the  oblique 
shock  i3  s'  1,  giving  large  expansion  of  the  nozzle  primary  flow. 

With  .  t  addition  there  are  several  changes  to  fne  flow  in  region  (2). 
iressure  ir  ncreased.  The  flow  is  deflected  to  larger  values  of  9.  The  Mach 
number  is  1  ered  in  region  (3)  due  to  heat  addition.  Either  a  weak  or  a  strong 
solution  of  he  oblique  shock  occurs.  When  the  strong  shock  solution  occurs, 
the  externa  flow  downstream  of  the  oblique  shock  is  high  subsonic.  Expansion 
waves  ref 1 .  ed  off  the  plug  surface  are  not  transmitted  into  the  external 
flow.  Hsnc  waves  are  trapped  and  are  reflected  as  compression  waves.  The 
compress irr  ives  maintain  a  high  plug  pressure.  Having  looked  at  the  inter¬ 
action  qua!  itively,  let  us  now  examine  the  results  of  a  sample  quantitative 
calculation. 

W-.ve  Tract  jr  ;  Quantitative  Example 


1 


The  sara  geometry  plug  nozzle  was  selected  as  that  tested  by  Harrington. 

.0 


(18) 


The  primary  *  ap  forming  the  boattail  was  turned  inward  at  an  angle  of  l7  .  A 
frees tream  V  h  number  of  1.31  was  chosen  for  the  calculations;  one  reason  for 
selecting  th  value  of  Mach  number  was  to  keep  the  external  flow  supersonic 
throughout  •  n  th>  was  no  heat  addition.  A  nozzle  pressure  ratio  of  6  was 
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chosen.  According  to  Figure  5,  at  maximum  thrust  a  higher  NPR  would  be  more 
appropriate  for  Mq  =  1.31.  NPR  of  6  is  characteristic  of  cruise.  In  order  that 
standard  gas  tables  could  be  used,  the  ratio  of  heat  capacities,  Y,  was  set 
equal  to  1.40. 

To  obtain  the  flow  field,  a  planar,  two-dimensional,  finite,  wave  calcula¬ 
tion  procedure  outlined  in  Section  12.9  of  Liepmann  and  Roshko^"^  was  used. 

(13} 

When  the  zone  of  influence,  as  defined  by  Oswatitsch  *  in  Section  3.25  of  his 
book,  spans  a  distance  small  compared  to  flow  radius,  the  two-dimensional  planar 
procedure  can  give  quite  accurate  answers  for  local  values  of  flow.  Near  the 
plug  nozzle  throat,  this  procedure  should  give  small  errors;  however,  near  the 
plug  tip,  sizeable  error  should  be  expected.  In  addition,  the  change  in  wave 
angle  at  the  intersection  of  finite  waves  was  neglected. 

Results  of  the  calculation  are  3hown  in  Figures  38  and  39.  Figure  38  shows 
the  wave  geometry.  The  dashed  line  is  the  slipstream  between  the  external  and 
the  nozzle  flows. 

Consider  the  Er  wave*  separating  regions  (l)  and  (2).  It  is  reflected 

as  an  expansion  wave  from  the  plug,  becoming  an  E^  wave.  The  Ej  wave  reflects 

from  the  slipstream  immediately  downstream  of  region  (4).  The  value  of  the 

reflection  coefficient  is  +  0.8.  This  means  the  reflected  wave  is  E  with  a 

r 

strength  of  0.  S  of  that  of  the  incident  Ej  wave.  At  the  next  reflection  from  the 
slipstream,  the  coefficient  is  *  0.12,  which  means  continued  expansion.  At  the 
next  reflection,  the  coefficient  is  -  0.08;  there  is  an  extremely  weak  com¬ 
pression  wave  reflected  from  the  slipstream.  As  a  consequence  of  the  expansion, 

*E  means  expansion,  and  C  means  s  compression  wave.  Subscript  r  is  for  a 
right  running  wave,  whereas  subscript  l  is  for  a  left  running  wave. 
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the  pressure  along  the  plug  drops  below  ambient  as  shown  in  Figure  39 .  Fressure 
below  ambient  on  the  plug  results  in  a  negative  plug  thrust.  Integration  of  the 
curve  shown  in  Figure  39  yields:  a  value  of  the  ratio,  plug  thrust  to  ideal  nozzle 
thrust,  equal  to  -  0.097.  Direct  comparison  of  this  calculation  with  Harrington's 
data  is  not  possible  since  all  his  tests  at  M  =  1.31  were  run  with  a  shroud 
extending  partially  over  the  plug.  However,  the  value  of  plug  thrust  has  the 
correct  algebraic  sign  and  the  correct  order  of  magnitude.  In  fact,  it  is 
better  than  the  correct  order  of  magnitude;  it  is  probably  within  50  per  cent 
of  being  correct. 

Now  let's  add  heat  to  the  region  near  the  boattail  as  shown  in  Figures  37 
and  40.  The  length  of  the  boattail  is  L.  The  thickness  of  the  heat  addition 
region  is  l.  The  pressure  increase  due  to  heat  addition  is  obtained  from 
Equation  (41).  Evaluation  of  the  integral  yields 


(136) 


Equation  (136)  can  be  solved  for  the  heat  addition  h,  BTU/volume  second,  and 
then  multiplied  by  Li. 


(13^) 


Eultj.plication  by  2*r,  i.e.,  the  circumference  of  the  boattail,  leads  to  the 
time  rate  at  which  heat  is  released.  Let  H  be  the  heating  value  of  the  fuel 
and  m^,  the  fuel  flow  rate.  The  heat  release  rate  is  given  by 

V  ■  Ms  <’ 


(138) 


The  specific  fuel  consumption,  SFC,  is  given  by 


3600m. 


SFC  = 


(139) 


FIGURE  40.  HEAT  ADDITION  REGION  AT  BOATTAIL. 
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when  AD  is  the  change  in  boattail  drag.  The  change  in  boattsil  drag  is 
AD  =  Lp’  sin  a;  the  angle  o  is  defined  in  Figure  40.  Combining  Equations 
(138)  and  (l79)  along  with  the  expression  for  AD  gives 

3600aS 


SFC  = 


(HO) 


(y  -  1)MH  sin  a 

Inserting  numerical  values,  a  =  1000  ft/sec,  p  =  -^M2  -  1  =  1.615,  Y  =  1.4, 

M  --  1.92,  H  =  20,000  BTU/lb^  and  t  =  17°,  and  dividing  by  778  ft-lbf/BTU 
leads  to  SFC  =  1.685  per  hour.  This  is  not  particularly  an  exciting  value  for 
SFC;  however,  for  the  strong  solution  the  heat  addition  causes  wave  trapping  and 
significantly  increases  plug  thrust.  Note  that  SFC  is  not  dependent  on  the 
amount  of  heat  added.  Less  heat  gives  less  thrust;  more  heat  gives  propor¬ 
tionately  more  thrust. 

Now  let's  look  at  the  flow  on  the  plug.  To  start  the  solution,  the  slip¬ 
stream  angle,  9,  must  be  determined.  Figure  41  is  a  plot  of  the  curves  showing 
pressure  in  the  external  flow  as  a  function  of  9  and  the  pressure  of  the  nozzle 
flow  also  as  a  function  of  9.  The  ordinate  is  pressure  behind  the  oblique 
shock  (or  downstream  of  the  expansion  fan  for  the  nozzle  flow)  divided  by 
freestream  stagnation  pressure.  There  are  two  possible  solutions;  the  strong 
solution  was  chosen  for  this  example.  Values  are  9  =  ♦  0.3°,  subsonic  external 
flow,  and  a  shock  wavs  angle  of  77°,  Knowing  9  at  the  boattail-slipstream 
junction  permits  one  to  start  the  solution  using  finite  waves. 

The  wave  geometry  along  the  plug  is  shown  in  Figure  42.  Since  the 
reflection  coefficients  at  the  slipstream  have  the  value  -  1.0,  the  flow 
becomes  periodic.  The  pressure  along  the  plug,  which  is  shown  in  Figure  43,  is 
everywhere  larger  than  ambient  pressure.  Integration  of  the  pressure  gives 
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If  heat  is  added  so  that  p’/p^g  =  .027  on  the  boattail,  the  value  of  SFC 
becomes  less  than  0.1  per  hour  In  Equation  (139),  AD  would  now  include 
both  the  decrease  in  boattail  drag  and  the  increase  in  plug  thrust. 

Wave  Trapping :  Discussion  of  Numerical  Example 

In  the  course  of  making  a  numerical  experiment,  there  are  numerous 
assumptions  and  decisions  to  make.  Having  completed  this  part  of  an  analysis 
of  the  influence  of  heat  addition,  one  has  new  perspective.  The  decrease  in 
boattail  drag  and  the  corresponding  SFC  are  fairly  straightforward  calcula¬ 
tions,  and  the  value  of  SFC  =  1.685  is  probably  fairly  accurate. 

Considering  the  plug  thrust  increment  due  to  heat  addition,  it  is  apparent 
that  a  crossroads  occurs  at  the  point  of  taking  a  weak  or  strong  solution. 

See,  once  again,  Figure  41.  For  this  example,  the  strong  solution  was  taken. 

The  strong  solution  gives  a  subsonic  external  flow  at  least  initially. 

That  fact,  of  course,  is  the  basis  for  wave  trapping.  Somewhere  along  the  plug, 
the  external  flow  may  accelerate  from  subsonic  flow  to  supersonic  flow.  Once 
tho  external  flow  becomes  supersonic,  the  waves  are  no  longer  trapped  in  the 
nozzle  flow.  Determination  of  the  subsonic  flow  and  the  change  to  supersonic 
flow  is  extremely  difficult. 

If  one  had  taken  the  weak  solution  instead  of  the  strong  solution,  the 
pressure  distribution  along  the  plug,  at  least  according  to  current  estimates, 
would  be  qualitatively  similar  to  Figure  39. 

Changing  NFR  raises  the  curve  labelled  "NOZZLE  FLOW"  in  Figure  41, 

As  the  nozzle-flow  curve  is  raised,  the  weak  and  strong  solutions  converge  to  a 
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single  solution  at  the  point  of  tangency.  See  Figure  44(a)  and  44(b), 

This  single  solution  is  a  strong  shock  solution.  Increasing  NFR  beyond  the 
point,  of  tangency  raises  the  curve  for  "NOZZLE  FLOW"  so  that  an  intersection 
of  the  nozzle  flow  curve  and  external  flow  curve  does  not  occur.  The  oblique 
shock  probably  moves  upstream  to  the  point  where  the  boattail  is  formed. 

This  case  is  illustrated  in  Figure  44(c).  This  would  give  high  pressure  on 
the  boattail  but  probably  would  give  low  plug  thrust. 

Adding  heat  on  the  boattail  has  two  influences;  (l)  Due  to  change  in  9, 
the  curve  for  external  flow  in  Figure  41.  is  shifted  to  the  right;  and  (2)  due 
to  increased  pressure  ahead  of  the  shock  wave,  the  curve  is  moved  vertically. 
This  is  illustrated  in  Figure  45  (a).  For  large  NFR  and  small  heat  addition, 
there  may  be  a  solution  as  shown  in  Figure  45(c).  Additional  heat  release 
causes  a  3hift  in  the  curve  making  possible  both  weak  and  strong  solutions. 

See  Figure  45(d). 

The  validity  of  the  flow  shown  in  Figures  42  and  4"!  is  doubtful  due  to 
the  fact  that  the  external  flow  probably  accelerates  from  subsonic  to  super¬ 
sonic.  Tne  strong  solution  may  not  occur.  Consequently,  one  is  skeptical 
of  the  SFC  =  0.1.  Additional  analysis  is  necessary  to  obtain  the  plug 
pressure  distribution  with  heat  addition.  Both  the  weak  solution  and  refine¬ 
ment  of  the  strong  solution  should  be  considered.  One  would  expect  tie  heat 
addition  to  have  a  favorable  influence  on  plug  thrust. 

Base  Flow  Froblem  with  Hnat  Addition:  Planar  Flow 

One  application  of  external  burning  which  has  been  investigated  both 
experimentally  and  theoretically  is  the  modification  of  base  pressure  by  heat 
addition.  The  heat  may  be  added  in  the  base  recirculation  zone,  in  the  viscous 
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9  DISTANCE  ALONG  PLUS  DISTANCE  ALONG  PLUG 

I  STRONG!  I  WEAK  i 


(a)  TWO  INTERSECTIONS-,  WEAK  AND  STRONG  SOLUTIONS. 


(STRONG  SAME  AS  ILLUSTRATED  ABOVE) 


(b)  TANGENT,  STRONG  SOLUTION. 


/ 


OBLIQUE  SHOCK 
RECIRCULATION  ZONE 


DISTANCE  ALONG  PLUG 


(c)  NO  INTERSECTION;  SHOCK  UPSTREAM  OF  BOATTAIL. 

FIGURE  44. 

INFLUENCE  OF  CHANGING  NOZZLE  PRESSURE  RATIO  WITH  FIXED  HEAT  ADDITION. 
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(a)  INFLUENCE  OF  INCREASING  HEAT  ADDITION  ON  EXTERNAL  FLOW. 


(b)  NO  INTERSECTION  WITH  SMALL  HEAT  ADDITION. 


PRESSURE  DISTRIBUTION 
AND  FLOW  SIMILAR  TO 
FIGURE  44(b) 


(c)  TANGENT  STRONG  SOLUTION. 


PRESSURE  DISTRIBUTION 
AND  FLOW  SIMILAR  TO 
FIGURE  44 (a) 


<d)  TWO  SOLUTIONS?  WEAK  AND  STRONG. 


FIGURE  45. 

INFLUENCE  OF  CHANGING  NEAT  ADDITION  WITH  FIXED  NOZZLE  PRESSURE  RATIO. 
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shear  layer,  or  a  the  inviscid  flow  adjacent  to  the  shear  layer.  In  1967, 

( lo ) 

Billig'’1  '  prepared  an  excellent  summary  of  work  to  date.  More  recently, 

1  v 

Heberts'  '  has  -ritten  a  survey  on  the  subject. 

It  is  wort; -  ;ile  to  examine  the  base  flow  from  an  integral  point  of  view. 

This  is  illustr  ->!d  in  Figure  46.  The  control  volume  for  the  momentum  theorem 

is  shown  as  das:  lines  with  various  .urfaces  identified  by  a,  b,  c,  .  .  .  . 

(15) 

Using  notation  of  Chapter  7  of  Liepraann  and  Roshko,  the  momentum  theorem  is 


J'pui(^nJdA  «■  Jpn^ 


dA  =  0 


(HD 


On**  car.  add  an  in  ■  egral  of  ambient  pressure  over  the  control  volume.  For 
base  pressure  i  =  L. 


ml(ulni  ♦  u2n2)dA  + 


(p  -  Pw)nxdA  =  0 


(142) 


For  the  control  v  .ume  shown  in  Figure  46,  there  are  four  integrals. 


r  f  P  f 

(Pb  -  PM)(“  1  dA  +  pu^(-  l)dA  +  I  (p  -  P^rijdA  + 


pu^dA  =  0  (143) 


be 


cd 


de 


(pb  -  pJdA  = 


nb 


pu^dA  ♦ 


de 


-  pu^dA  ♦  ^  (p  -  p^^dA  (144) 

0  bc  ©  cd  © 

Heai  addition  in  t. e  inviscid  flow  above  the  edge  of  the  shear  layer  will  alter 

the  shape  of  tne  0  reamline  from  c  to  d.  Changing  streamline  cd  will  not 

influence  integral  '2^  above.  It  is  not  readily  apparent  how  integral  (T) 


1CT7 


S^SSfw^S??? 


changes  due  to  modification  of  streamline  cd.  Of  course  integral  ^3^  is 


dire^ly  influenced. 


To  gain  some  insight  into  the  anticipated  change  in  integral  one 


can  use  the  continuity  equation  which  is 


pux 

*w 


•u^dA  «  1  pu^dA 


(145) 


Since  integral  over  surface  be  doesn't  change  when  streamline  cd  is  changed, 

this  means  the  integral  over  de  does  not  change.  Any  variation  in  integral 

2 

of  Equation  (144)  is  due  to  variations  in  p  and  u^;  however,  Equation  (145) 
tells  us  any  change  in  p  is  balanced  by  a  compensating  change  in  u^» 

If  the  flow  above  the  edge  of  the  shear  layer  illustrated  in  Figure  46 


is  inviscid  and  if  only  one  family  of  characteristics  is  significant,  then 


one  can  relate  local  streamline  slope  to  pressure.  To  demonstrate  the  influence 


of  the  shape  of  streamline  cd  on  base  pressure,  a  sample  problem  was  worked. 


This  is  illustrated  by  Figure  47.  Streamline  cd  is  used  with  the  same  lateral 


displacement  but  with  large  turning  angle  in  the  upper  diagram.  The  pressure 


distribution  along  the  streamline  is  shown  in  the  graphs  below  each  streamline. 


From  the  curves  one  can  evaluate  the  integral  or  equivalent  series 


f  p  V""'  Pj 

=  1  — r»,dA  =  }  -*AA, 

J  pT  1  Z_y  pT  > 


(146) 


where  AA^  is  the  element  of  area  projected  in  the  direction.  Numerical 
values  for  I  are  12.8  for  upper  curve  and  21.7  for  lower  curve.  Relating 


Equation  (146)  to  Equation  (144)  gives 


— n.sii 


SIAM  AUlWMMl  (MtM»  oeKNMMM  uttMuAWl. 


(-  -  “)dA  = 


P<p 


(p  ~  P*)nidA  =  P-p^1  “  52.8) 


(147) 


cd 


cd 


The  sharp  turning  curve  gives  a  value  of  -  40p^  for  integral  of 


Equation  (144)»  whereas  the  gradual  turning  curve  gives  -  31. Ip^.  The  base 


pressure  would  be  greater  with  the  gradual  turning  streamline  cd.  It  is 
obvious  from  this  analysis  that  one  wants  to  maintain  a  high  pressure  on 
streamline  cd  to  increase  base  pressure.  The  heat  addition  zone  shown  in 
Figure  48  will  do  that. 

Now  let's  return  to  the  problem  of  how  much  pressure  increase  results 
from  the  heat  addition  zone.  Consider  a  heat  addition  zone  located  above  a 
slipstream  as  shown  in  Figure  49.  This  is  a  model  to  represent  the  flow 
depicted  in  Figure  48.  One  replaces  the  lateral  gradient  in  Mach  number 
with  a  slipstream  having  supersonic  flow  above  and  subsonic  flow  below.  V?ith 
a  lateral  Mach  number  gradient,  as  in  Figure  48,  the  characteristics  reflect 
from  the  sonic  line  in  the  form  of  a  cusp. 

A  paragraph  will  now  be  used  to  explain  the  notation  E^,  C E^,  C^. 

The  symbols  E  and  C  are  for  Expansion  and  Compression  waves.  The  subscripts 
r  and  l  represent  right  running  and  left  running  waves  respectively. 

In  Figure  49 the  heat  addition  zone  has  a  leading  edge  swept  at  angle  jj. 
This  gives  a  finite  wave  Cr.  The  pressure  along  the  slipstream  remains  at 
pressure  p^.  The  waves  from  the  heat  addition  zone  and  the  reflected  waves 
turn  the  flow  while  keeping  pressure  constant.  One  way  to  describe  this  is  to 
say  that  the  compression  waves  from  the  heat  addition  zone  have  cancelled  the 
expansion  waves  which  would  have  originated  due  to  kink  in  streamline  at 
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FIGURE  48. 

INTERACTION  OF  HEAT  ADDITION  ZONES  WITH  BASE  FLOW. 


point  0.  The  turning  angle  is  2©.  at  point  (0  or  twice  the  turning  angle  due 
to  C^.  The  reflection  has  effectively  doubled  the  influence  of  the  waves  from 
the  heat  addition  zone. 

The  doubling  of  effectiveness  can  be  seen  in  the  following  analysis:  In 
free  unbounded  space  the  turning  angle  is  giver,  by 


tan  0=0  = 


U  ♦  u'  ~  U 


(148) 


Equation  (148)  is  obtained  from  geometry  and  Equation  (40).  From  Equation  (41 ) 


(149) 


where  symbol  b  is  defined  by  Figure  49.  Due  to  turning  of  flow  by  an 
amount  29,  the  increment  of  pressure  decrease  would  be 

pH  -  _  &Mjb8 


(150) 


From  the  flow  field  analysis  illustrated  in  Figure  49,  the  following 


equation  must  be  true 


2p*  ♦  prt  =  0 


(151) 


Combining  Equations  (148),  (149),  (150),  and  (151),  one  can  verify  that 
Equation  (151 )  is  in  fact  true. 

The  mechanism  for  altering  base  pressure  by  heat  addition  in  the  adjacent 
inviscid  flow  is  apparent  if  one  retraces  the  discussion  associated  with 
Figures  46,  4?,  and  48.  To  summarize,  the  streamline  cd  must  turn  inward 
toward  the  axis  or  plane  of  symmetry.  Without  heat  addition,  this  turning 
causes  the  pressure  to  drop.  Sharp  turning  gives  low  pressure  and  low  base 
pressure;  see  Equation  (144)  and  Figure  47*  Heat  addition  allows  this  inward 


turning  without  loss  of  pressure  as  shown  in  Figures  48  and  49.  Strahle’s 
analysis  suggests  a  new  way  to  overcome  base  drag. 
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Base  Flow  Problem  with  Heat  Addition ;  Axisymmetric 


(22)  (2^) 

The  analysis  of  Fein'  '  and  Strahle'  '  considered  a  two-dimensional 


planar  geometry,  i.e.,  wing-like  objects.  Projectiles  are,  of  course, 
axisymmetric.  This  section  examines  some  of  the  changes  to  be  expected  when 


the  planar  results  are  replaced  by  calculations  for  the  axisyiriietric  case. 

.(15) 


The  method  outlined  in  Liepmann  and  Roshko'  '  was  used  to  calculate  the 
flow  field  for  axisymmetric  heat  addition.  To  verify  the  procedure  for 
handling  the  characteristics  near  the  axis,  a  model  of  a  spherical,  radially 
expanding  flow  was  used.  The  results  of  that  calculation  are  shown  in  Figure  5CI 
Since  the  flow  is  known,  one  knows  Mach  number  at  any  radius  R.  The  Mach 
number  calculated  from  the  characteristics  solution  is  shown  in  the  third 
column  of  the  data  summary  in  Figure  5&  The  Mach  number  based  on  the  value 
of  R  and  flow  area  is  shown  in  the  first  column. 


A  planar  radial  flow  has  an  area  variation  linear  in  R,  whereas  the 

,2 


axisymmetric  radial  flow  has  a  R  dependence  for  area.  This  can  be  represented 
by 

A  A  A T3  A  A  ATI 

(152) 


AA  _  AR 
A  "  R 
Planar 


AA  _  ~AR 
Axisymmetric 


Partial  solutions  of  the  flow  field  with  heat  addition  were  obtained  for 
two  cases.  Both  cases  had  an  annular  heat  addition  zone  as  illustrated  in  the 
lower  left-hand  corner  of  Figure  51.  The  heat  addition  was  sufficient  to  cause 


a  maximum  turning  angle  of  9  .  =  -  6  at  a  radius  of  6  inches.  In  the  first  case, 

max 


the  heat  zone  was  bounded  by  6  ^  r  4  ?  and  0  ^  z  ^  3.  In  the  second  case,  it 
was  6  ^r  ^8  and  0  ^  z  ^  T.  For  the  first  case,  shown  in  Figure  51,  Mq  *  V* 
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FIGURE  50.  CHARACTERISTICS  SOLUTION  FOR  AXISYMMETRIC 
RADIAL  EXPANSION. 


For  the  second  case,  Mq  was  chosen  as  2,  and  the  radial  extent  of  the  heat 
addition  zone  was  increased  to  spread  out  the  waves.  This  was  done  to  avoid 
merging  of  the  characteristics  to  form  a  shock  wave. 

Pressure  versus  radius  is  plotted  in  Figure  51*  The  mid-ehara^eristic, 
which  has  points  1,  2,  and  3,  has  an  imperceptible  increase  in  pressure  as 
radius  decreases.  The  characteristics  converge  to  form  a  shock  wave.  When  this 
happens,  the  techniques  outlined  in  the  bookv<  and  NACA  TIT  '  by  Ferri  are 
appl ied . 

Along  the  characteristic,  which  starts  at  r  =  6"  with  6  =  -  6°,  the 

pressure  rises  slightly  to  a  radius  of  r  =  2,f.  Between  0  ^  r  ^  2,  the  pressure 
increases  rapidly.  Due  to  the  large  mesh  size,  the  pressure  variation  at 
points  12,  13,  and  14  may  not  be  correct.  Certainly  the  end  point,  number  14, 
is  correct.  The  dashed  line  represents  the  curve  that  would  probably  • 
obtained  with  finer  mesh  size  and  a  more  accurate  treatment  of  the  shock  waves. 
Note  that  for  =  2,  even  a  normal  shock  wave  is  nearly  isentropdc,  which 
permits  treatment  of  shock  waves  as  finite  isentropic  compression  waves. 

The  conclusions  from  Figure  51  are  as  follows: 

(a)  In  the  axisymmetric  geometry,  compression  waves  tend  to  merge  more 
quickly. 

(b)  The  increase  of  pressure  with  decreasing  radius  is  very  small  until 
the  axis  is  approached. 

Now  let's  look  at  the  second  example  which  was  chosen  to  avoid  a  shock  wave. 

The  characteristics  for  the  second  case  are  shown  in  Figure  52.  Increasing 
the  spacing  of  the  waves  did  not  completely  avoid  a  shock  w«ve.  A  normal  shock 
extends  from  points  13  are!  14  to  the  axis. 


- x — f—  FIGURE  52.  INITIAL  CHARACTERISTICS  IN  REGION 

HEAT  ADDITION  ZONE  T  1  INTERNAL  TO  ANNULAR  HEAT  ADDITION. 


The  pressure  as  a  function  of  radius  is  shown  in  Figure  53  for  the  raid-wave 
(wave  starting  at  r  =  6  with  0  =  -  3°)  and  the  terminal  wave  (wave  starting  at 
r  =  6  with  ©  =  -  6°).  As  before  in  the  previous  example,  the  pressure  does  not 
increase  significantly  until  very  close  tc  the  axis.  This  is  true  for  both  the 
mid-wave  and  the  terminal  wave.  The  pressure  is  snown  as  a  dash-dot-dash  line 
between  point  14  and  the  axis  since  this  was  not  calculated  in  detail. 

Figure  54  gives  information  about  the  streamlines  and  flow  deflection. 

Right  on  the  axis,  of  course,  the  streamline  is  the  axis.  There  is  large 
tur*  i  near  the  axis  giving  small  radius  of  curvature  for  the  streamline. 

This  balances  the  large  radial  pressure  gradient  near  the  axis. 

The  conclusions  from  the  second  case  are  the  same  as  the  first  ca3e. 
Preliminary  Comments  Based  on  a  Cursory  Look  at  Axisymmetric  Case 


The  change  in  base  drag  is  given  by 

ADb  =  ApbAb 


(153) 


Each  of  the  two  terms,  Apb  and  Ab»  will  now  be  discussea.  The  results  shown 
in  Figures  51,  52,  and  53  suggest  that  Apb  may  not  be  Increased  significantly 
as  a  direct  result  of  changing  geometry  from  planar  to  axisymmetric.  However, 
the  results  shown  in  Figures  51,  52,  and  54  more  closely  apply  to  pressure  at 
the  dividing  streamline  and  not  to  base  pressure.  The  flow  in  the  base  region 
may  magnify  the  pressure,  just  as  the  initial  pressure  due  to  heat  addition  at 
r  =  b  .s  increased  slightly  due  to  axisymmetric  geometry. 

Compression  originating  with  the  heat  addition  will  be  reflected  at  the 
sonic  line.  See  Figure  48  and  point  of  Figure  55.  These  compression  waves 
will  h**  reflected  from  the  wake  of  the  heat  addition  zone;  see  point  (b)  of 
Figure  55.  The  reflection  may  be  negative,  i. e. ,  changing  the  compression  waves 

to  expansion  waves.  This  negative  reflection  would  not  be  favorable. 
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Now  let's  look  at  the  area  terra.  Consider  the  seme  amount  of  heat,  Q, 
released  in  both  planar  and  axisymmetric  case.  In  the  axisymmetric  case,  Q  is 
distributed  along  a  circle  of  radius  r ^  giving  k/2^r^,  BTU/length  3ec.  In  the 
planar  case,  Q  is  distributed  along  two  lines  of  length  mr^.  One  line  i3  at  the 
top,  and  the  other  at  the  bottom.  Sea  Figure  56.  For  the  planar  case,  the 
heat  release  is  Q/^r^,  BTU/length  sec,  which  is  identical  to  the  axisymmetric 
case.  The  area  influenced  by  the  heat  release  is  in  the  ratio 

2 

base  area;  axisymmetric  case  mr.  r. 

-  -  _  =  ■!»—  /it  l  \ 

base  area;  planar  case  ^rlT2  ^r2  v 

Making  r^/r^  small  helps  increase  pressure  at  r^  due  to  heat  addition;  see 
Figure  56.  Making  r-j/r^  lar6e  makes  the  base  area  very  small  in  axisymmetric 
case  as  compared  to  the  planar  case.  One  factor  helps;  the  other  hurts. 

As  a  final  remark,  it  should  be  stated  again  that  these  are  preliminary 
thoughts  and  conclusions. 

Spin  Recovery  Wing  Heat  Addition 

Some  modern  fighter  and  attack  aircraft  can  be  flown  into  spins  of 
such  nature  that  special  equipment  is  needed  for  recovery.  The  drag 
parachute  may  be  deployed  to  lift  the  aircraft  tail.  Before  normal  flight 
is  restored,  the  parachute  must  be  jettisoned.  An  alternate  method  of 
obtaining  a  torque  on  the  aircraft  would  be  heat  addition.  Figure  57  illustrates 
a  modern  fighter  in  a  spin.  The  engines  run  normally  and  produce  thrust.  By 
spraying  fuel  into  the  exhaust  stream,  the  heat  addition  deflects  the  jet  down~ 
ward.  A  torque  is  developed  which  tends  to  lift  the  tail.  Let's  compare  a 
parachute  recovery  with  heat  addition. 
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Assume  the  aircraft  has  a  moment  of  inertia  of  6  x  10  lb  ft  about  the 
pitch  axis.  Assume  the  installed  thrust  is  35,000  lb  with  afterburner.  If 
the  flow  can  be  deflected  by  6°  with  heat  addition,  the  force  tending  to  lift 
the  tail  is  4000  lb.  If  the  aircraft  has  a  moment  arm  of  22  feet,  the 
restoring  torque  is  8.8  x  lO^lb  ft.  Time  required  to  point  the  aircraft 
downward,  i.e. ,  turn  90°  in  pitch,  is 

t  =  3  46  sec 

If  the  sink  rate  is  140  ft/sec,  the  recovery  altitude  is  6400  ft. 

One  can  estimate  the  amount  of  heat  to  be  added  for  a  deflection  of  6° 
using  Equations  (119)  through  (129).  The  temperature  of  exhaust  must  be 
increased  hy  30  per  cent. 

Now  consider  a  parachute  of  10  ft  diameter  with  a  drag  coefficient  of  2. 

The  drag  at  a  sink  rate  of  140  ft/sec  will  be 

D  =  -rpt^AGjj  =  4200  lb 

This  is  comparable  to  the  force  due  to  thrust  deflection.  Performance  of 
parachute  and  heat  addition  is  comparable. 

SUMMARY  AND  CONCLUDING  DISCUSSION 

Use  of  an  area  rule  or  quasi  area  rule  implies  optimization.  The  area 
rule  can  be  used  to  translate  an  optimum  axisyrametric  body,  e.g. ,  von  Karman 
ogive,  to  a  three-dimensional  body.  The  optimum  axisymmetric  body  shape  is  found 
by  application  of  the  calculus  of  variations  to  the  source  distribution  function, 
f(x).  When  there  is  heat  addition,  there  is  a  new  function  to  optimize,  Q(x,r); 
see  Figure  55.  Both  f(x)  and  Q(x,r)  must  be  simultaneously  optimized.  The 


axisyrametric,  small  perturbation,  solution  for  heat  addition  has  never  been 
obtained;  it  is  an  essential  building  block  of  the  quasi  area  rule. 

When  the  body  develops  lift  with  aerodynamic  circulation,  there  is  another 
complexity  to  the  analysis.  Vorticity  due  to  lift  must  be  included  in  the  analysis. 
There  may  be  an  optimum  balance  between  lift  generated  aerodynamically  and  lift 
generated  by  heat  addition. 

To  apply  the  optimization  technique  to  turbojets,  ramjets,  or  other 
internally  burning,  air  breathing,  engines,  a  different  analytical  technique 
is  required.  A  combined  thrust-drag  technique  was  briefly  introduced.  This 
method  uses  an  energy  disc. 

By  discussing  both  one-dimensional  and  two-dimensional,  planar,  heat 
addition,  the  connection  between  heat  fronts,  combustion  fans,  and  one-dimensional 
heat  addition  was  demonstrated.  Heat  addition  in  the  planar  case  may  be  represented 
either  by  line  sources  or  by  a  volume  distribution  function,  Q(x,r).  In  several 
ways  heat  addition  is  similar  to  the  flow  generated  by  a  solid  body.  The 
pressure  cn  the  flat  plate  of  Figure  32  is  the  same  as  that  caused  by  a  wedge 
of  angle  0.  Fressure  and  velocity  perturbations  are  identical  for  both  solid 
body  and  heat  addition.  Differences  exist  between  flow  over  solid  bodies 
and  flow  through  a  region  of  heat  addition.  One  difference  is  drag. 

There  are  three  approaches  to  calculate  thrust  minus  drag  in  an  inviscid 
flow:  integration  of  pressure  over  the  body  surface,  momentum  control  vol:une, 
and  wave  drag.  Integration  of  pressure  is  the  most  direct  of  the  methods. 

Figures  22  and  27  illustrate  the  wave  energy  and  momentum  approaches  for  a 
thrusting,  lifting,  planar  airfoil.  Figures  20,  23,  and  26  illustrate  the  wave 
and  momentum  approaches  for  an  isolated  heat  addition  zone  in  an  infinite  medium. 

The  merit  of  a  particular  method  depends  on  the  geometry. 

128 


**&Z5Wi 


Applications  of  external  burning  include  lifting,  thrusting  airfoils,  base 
or  boattail  pressure  modification  for  an  exhaust,  development  of  base  thrust 
for  projectiles,  and  spin  recovery.  External  burning  for  airfoils  and  projectiles 
has  been  studied  extensively  both  theoretically  and  experimentally.  A  new 
phenomenon  occurs  for  the  exhaust  problem;  this  is  wave  trapping.  A 

i 

preliminary  analysis  indicates  heat  addition  for  spin  recovery  is  comparable 
to  use  of  a  parachute. 
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